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ABSTRACT 

The  problem  of  e  partial  or  aupercavitating  hydrofoil  in  linear 
theory  ia  formulated  in  terms  of  integral  equationa  of  unknown  vortex 
and  cavity  source  distributions.  The  general  problem  ia  decomposed 
into  one  of  camber,  one  of  thickness,  and  one  of  angle  of  attack. 

The  general  solution  is  given  in  terms  of  integrals  of  known 
functions  over  the  cavity  length  for  partial  or  over  the  chord  length 
for  aupercavitating  hydrofoils.  The  numerical  scheme  to  compute  the 
integrals  i6  shown  to  be  very  accurate  and  insensitive  to  the  variables 
of  the  problem. 

Cavity  shapes  obtained  from  linear  and  non-linear  theory  are 
compared  for  certain  special  cases.  Linear  theory  ia  shown  to  be 
accurate  for  aupercavitating  sharp-nosed  hydrofoils  at  moderate  angles 
of  attack. 

The  effect  of  the  leading  edge  radius  on  the  partial  cavitation  of 
a  hydrofoil  is  introduced  by  incorporating  Lighthill's  correction  in 
the  linearized  formulation  of  the  problem.  The  results  for  certain 
special  cases  show  the  significant  role  of  the  leading  edge  radius  on 
the  cavitation  of  hydrofoils,  something  also  confirmed  from  experiments 
conducted  at  the  MIT  water  tunnel. 
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Figure  2.  Partially  cavitating  hydrofoil  -  Linear  theory 


density  of  the  fluid,  ve  define  the  cavitation  number  o  as: 


o 


2 


(2.1) 


The  cavity  thickness  h(x)  has  to  be  determined  as  a  part  of  the 
solution  to  the  problem.  The  point  on  the  hydrofoil  where  the  cavity 
starts  can  be  determined  by  the  Villat-Brillouin  condition  [4,40] 
according  to  which  the  curvature  of  the  cavity  and  foil  should  be 
continuous  at  that  point.  The  end  region  of  the  cavity  can  be 
sianilated  by  using  sny  cavity  termination  model,  for  example,  the 
modified  Kiabouchinsky  model  or  the  reentrant  jet  model.  A  review 
of  the  proposed  models  is  included  in  [35] . 

Linear  theory  can  be  applied  if  the  slopes  of  the  hydrofoil  and 

cavity  surfaces  with  respect  to  the  main  stream  are  of  the  order  of  our 

small  parameter  e ,  which  can  be  for  example  the  maximum  thickness 

over  chord  ratio.  In  that  case  the  boundary  conditions  are  applied  on 

the  slit  from  x-0  to  x*l  as  we  can  see  in  Figure  2.  If  n^(x)  and 

t^(x)  are  the  upper  and  lower  surfaces  of  the  hydrofoil  respectively, 

and  (u  ,  v  )  the  perturbation  velocity  flow  field,  then  the 
c  c 

linearised  boundary  conditions,  as  demonstrated  in  [34]  f  vill  be: 

The  kinematic  boundary  condition 

vc(x,0+)  -  ;  i<x<l  (2.2) 


vc(x,0-)  -  »  0  <x<1  (2.3) 

^Throughout  our  analysis  it  has  been  assumed  that  the  cavity  starts  at 

the  leading  edge  of  the  hydrofoil.  This  is  further  discussed  in 
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Section  22. 


22 


The  streetwise  perturbation  velocity  u^  on  the  slit  0<x<l  will  be 
given  as: 


u  (x,0+)  = 
w 


_  Yu(x) 

-  =  -v- 


(1.4) 


The  thickness  problem: 


The  thickness  distribution  t(x)  is  represented  by  the  source 
distribution  q^(x).  The  linearized  kinematic  boundary  condition  will 


U  .-L-!1  =  vT(x,0+)  = 
°°  2  3x  w 


T  qw(x) 

-  0-)  =  -~2 


(1.5) 


Therefore  the  source  distribution  will  be: 


q  (x)  =  U  -I1 
w  «  3x 


(1.6) 


and  the  streamwise  perturbation  velocity  u^  on  the  slit  0<x<l : 


uT(x,0+)  =  uT(x,0-) 
w  w 


_1  rl 

2tt  ’  t  r-v 


(1.7) 


2.  Partially  cavitating  hydrofoil  -  Linear  theory. 

Consider  the  same  hydrofoil  as  in  the  previous  section  but  now  with 
a  cavity  of  length  i  on  its  upper  surface  as  shown  in  Figure  2. 

If  pc  is  the  corresponding  pressure  inside  the  cavity  and  P  the 


„-v  %  s  % 


Figure  1.  Hydrofoil  in  non-cavitat 


PART  A.  PARTIALLY  CAVITATING  HYDROFOIL-LINEAR  THEORY 


1 .  Hydrofoil  in  non-cavitating  flow  -  Linear  theory. 

Consider  a  hydrofoil  of  chord-length  c“l ,  camber  distribution 
n  (x)  (including  the  angle  of  attack  °0  and  thickness  distribution 
t(x)  in  a  uniform  flov  stream  U^,  and  ambient  pressure  ,  as 
illustrated  in  Figure  1.  According  to  linear  theory  described  in  more 
detail  in  [28]  the  perturbation  velocity  flow  field  (uw>  v^) 
can  be  generated  by  a  distribution  of  vortices  Y ^(x)  and 
sources  q^(x)  on  the  slit  from  x*0  to  x*l.  Furthermore  the  general 
hydrofoil  problem  can  be  decomposed  in  the  following  problems: 

The  camber  problem: 

The  camber  distribution  n (x)  is  represented  by  the  vortex 
distribution  yw(x).  The  linearized  kinematic  boundary  condition 
will  be: 


U  ^  *  vC 

U°°*  9x  w 


(x,±0)  = 


1Y  U).dc,* 


w 


e-* 


0<x<l 


(1.1) 


The  Kutta-condition  at  the  trailing  edge  of  the  hydrofoil  will  be: 


Y  (Definite  (1.2) 

w 

Inversion  of  the  singular  integral  equation  of  Cauchy  type  (1.1)  in 
combination  with  (1.2)  will  give  us  the  unique  solution: 


Y„U> 


(1.3) 


*The  symbol  j  denotes  the  Cauchy  principal  value  of  the  singular 
integral. 
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super cavit at ing  hydrofoils  is  formulated  in  a  uniform  way  and  the 
solution  is  given  in  terms  of  integrals  of  known  functions.  Then  the 
non-linear  foil  thickness  effects  are  introduced  by  incorporating 
Lighthill's  correction  in  the  linearized  version  of  the  cavitating 
hydrofoil  problem.  That  leads  to  an  equivalent  linear  problem  with  the 
described  stresmwise  perturbation  velocity  varying,  rather  than 
constant,  over  the  cavity  length.  The  present  theory  is  found  to 
predict  a  higher  reduction  in  the  cavity  size  than  that  by  Uhlman,  by 
increasing  the  foil  thickness,  but  it  is  not  as  high  as  that  predicted 
by  Tulin  and  Hsu.  Finally,  the  present  theory  is  expected  to  be  easily 
implemented  into  the  cavitating  propeller  program  via  a 
three-dimensional  version  of  Lighthill's  correction.  This  is  expected 
to  result  in  a  negligible  change  in  the  computing  time. 


numerical  methods  have  also  been  developed  for  the  analysis  of  the 
cavitating  hydrofoil  or  ving  problem  in  linear  theory.  First,  the 
discrete  vortex  and  source  method  has  been  applied  to  a  cavitating  flat 
plate  [14].  Then  an  unsteady  lifting  surface  theory  has  been  produced 
for  a  supercavitating  ving,  by  using  a  vortex-source  lattice  technique 
[17].  Finally,  this  numerical  method  has  been  extended  to  predict  the 
unsteady  sheet  cavitation  of  marine  propellers  operating  in  non-uniform 
wakes  [24,19,21].  It  has  been  observed  that  the  numerical  linear 
theory  overpredicts  the  extent  and  volume  of  the  cavities  for  certain 
propeller  applications.  This  has  been  attributed  to  the  inability  of 
the  linear  theory  to  include  non-linear  foil  thickness  effects  on  the 
cavitation  of  hydrofoils. 

Tulin  and  Hsu  developed  their  short  cavity  theory  [36]  by 
considering  the  cavitating  flow  as  a  small  perturbation  on  the 
non-linear  fully  wetted  flow.  Thus  they  were  able  to  include 
non-linear  foil  thickness  effects  in  their  formulation.  Namely,  they 
predicted  that  by  increasing  the  thickness  of  a  partially  cavitating 
hydrofoil,  the  size  of  the  cavity  was  reduced  substantially  for  fixed 
flow  conditions. 

Uhlman  [38]  formulated  a  non-linear  numerical  method  for  the 
analysis  of  partially  cavitating  hydrofoils  by  using  surface 
singularity  techniques  and  by  applying  the  exact  boundary  conditions  on 
the  cavity  and  on  the  foil.  His  method  also  predicted  a  reduction  in 
the  size  of  the  cavity  as  the  foil  thickness  increased,  but  cot  as  high 
as  that  by  Tulin  and  Hsu.  However,  the  computing  time  of  his  numerical 
scheme  was  long  and  therefore  it  was  difficult  to  be  implemented 
efficiently  into  the  cavitating  propeller  program  [21]. 

In  this  thesis,  first  the  linear  theory  for  partial  and 
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of  using  the  hodograph  technique  to  deal  with  curved  body  boundaries, 
very  few  cases  have  been  treated  analytically  by  using  the  exact 
non-linear  theory. 

Linear  theory  vas  first  applied  by  Tulin  [34]  to  the  problem  of  a 
supercavitating  flat  plate  at  snail  angles  and  arbitrary  cavitation 
nunbers.  It  was  subsequently  extended  by  T.  T.  Wu  [42],  Geurst  [11], 
Parkin  [30],  and  Fabula  [8]  for  a  supercavitating  hydrofoil  of  general 
shape.  These  authors  have  worked  by  means  of  the  complex  velocity 
function,  and  have  given  expressions  for  the  cavitation  number,  the 
hydrodynamic  coefficients  and  the  cavity  volume,  in  terms  of  integrals 
of  known  quantities. 

The  partially  cavitating  hydrofoil  problem  has  first  been  addressed 
in  linear  theory  by  Acosta  [1],  Geurst  and  Timman  [9]  independently, 
who  gave  analytical  results  for  a  flat  plate  inclined  to  a  uniform 
inflow.  Geurst  [10]  also  formulated  the  linear  theory  for  partially 
cavitating  camber  meanlines  and  gave  analytical  results  for  a  parabolic 
camber  distribution  [12].  The  foil  thickness  was  first  included  by  Wade 
[41]  who  analyzed  the  partially  cavitating  planoconvex  hydrofoil 
problem. 

Banaoka  [16]  stated  the  linear  partial  and  supercavitating 
hydrofoil  problem  in  terms  of  integral  equations  of  unknown 
distributions  of  sources  and  dipoles  on  the  foil  and  on  the  cavity. 
Inversion  of  these  integral  equations  produced  integral  representations 
for  the  cavitation  number,  the  hydrodynamic  coefficients  and  the  slope 
of  the  cavity  surface  in  terms  of  the  shape  of  the  hydrofoil.  Be  also 
gave  series  representations  for  his  results  when  the  hydrofoil  shape 
could  be  expressed  in  terms  of  polynomials  of  the  chordwise 


coordinate 


INTRODUCTION 


Cavitation  is  the  phenomenon  in  which  a  cavity  of  gas  exists 
inside  a  flow  stream,  due  to  the  development  of  pressures  lover  than 
the  saturation  pressure  at  the  local  temperature. 

This  is  very  likely  to  occur  on  the  sunction  side  of  a  hydrofoil 
or  on  the  backface  of  propeller  blades.  Since  in  our  time  ships  become 
larger  and  propeller  loadings  bigger,  cavitation  becomes  more  and  more 
inevitable  without  any  considerable  sacrifice  in  propulsion  efficiency 
and  therefore  fuel  consumption. 

On  the  other  hand,  unsteady  propeller  blade  cavitation,  which  is 
the  most  usual  case,  is  undesirable.  One  of  the  reasons  is  the  chance 
of  failure  at  the  propeller  trailing  edge,  due  to  the  impulsive  loads 
during  the  collapse  stage  of  the  unsteady  cavity.  However,  the  most 
important  reason  is  that  unsteady  cavitation  is  the  main  contributor  to 
unsteady  hull  forces,  which  often  cause  severe  vibrations  and 
structural  failure  in  some  ships. 

In  order  to  design  an  efficient  propeller  with  a  low  level  of 
cavitation,  a  thorough  understanding  of  the  phenomenon  of  cavitation  is 
desirable. 

Cavitating  or  free-streamline  flows  have  been  studied  extensively 
during  the  last  century  with  increasing  interest.  The  first  problems 
to  be  addressed  were  those  of  flows  around  bodies  at  zero  cavitation 
number.  The  analysis  of  these  problems  has  been  achieved  via  the 
hodograph  technique  [3].  The  extension  to  non-zero  cavitation  number 
problems  has  created  a  lot  of  diversity  in  the  cavity  termination 
models.  The  most  prevailing  models  seem  to  have  been  the  Riabouchinsky 
model  [32]  and  the  re-entrant  jet  model  [6].  Due  to  the  difficulty 
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Leading  edge  radius  of  the  hydrofoil 


Cavitation  number  defined  as: 


P»-Pc 
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Inception  cavitation  number 
Cavitation  number  in  pure  linear  theory 
Thickness  distribution  of  the  hydrofoil  (see  Figure  1) 
Maximum  thickness/chord  ratio  for  the  hydrofoil 
Thickness  at  the  trailing  edge/chord  ratio  for  the 
hydrofoil 


Transformed  £  variable  defined  as  £  *  fc.  sin2(-|);  ()<£<«, 


and  0<<t><iT 

Discrete  4*  defined  in  (8.7) 
Function  defined  in  (6.6) 
Dummy  z  variable 


Superscripts 

c:  Refers  to  the  camber  problem 

t:  Refers  to  the  thickness  problem 

a:  Refers  to  the  flat  plate  problem  at  angle  of  attack 

+:  Refers  to  the  upper  hydrofoil  surface 

- :  Refers  to  the  lover  hydrofoil  surface 

1:  Corresponds  to  No.  1  hydrofoil  (see  Section  24) 

2:  Corresponds  to  No.  2  tested  hydrofoil  (see  Section  24) 

Corresponds  to  No.  3  tested  hydrofoil  (see  Section  24) 
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Vortex  distribution  for  the  non-cavitating  hydrofoil 
problem  (see  Figure  1) 

Maximum  of  the  vortex  distribution  for  NACA  a  ■  m  camber 
meanline  at  a  ■  (see  Figure  5a) 

Thickness  at  the  end  of  the  cavity  (see  Figure  32a) 

Increment  in  the  angle  of  attack 
Change  in  o  due  to  open  cavity  model 

Perturbation  parameter  for  the  cavitating  hydrofoil  problem. 

Also  ratio  of  the  minor  to  the  major  axis  for  an  ellipse 
Dummy  z  variable 

Camber  distribution  of  the  hydrofoil  without  including  the 
angle  of  attack  (see  Figure  1) 

Maximum  camber /chord  ratio  of  the  hydrofoil 
Camber  distribution  of  the  hydrofoil  including  the  angle  of 
attack  (see  Figure  1) 

Cavity  surface  distribution 

Lower  hydrofoil  surface  distribution  (see  Figure  2) 

Upper  hydrofoil  surface  distribution  (see  Figure  2) 

Transformed  n  variable  defined  from  the  expression  u  *  s.sin 

2  0 

for  partially  cavitating  or  n  ■  t.sin  (y)  for 

supercavitating  hydrofoils;  CK0<ti 

*  1 

Hydrofoil  slope,  defined  as:  e  (x)  - 

0  oX  • 

The  maximum  camber/chord  ratio  for  a  parabolic  camber 
zieanline,  doubled  (see  Figure  4b) 

Maximum  thickness/chord  ratio  for  a  biconvex  hydrofoil  (see 
Figure  4c) 


o: 


Dummy  x  variable 
Density  of  the  fluid 


x-conponent  of  Che  perturbation  velocity  with  respect  to  U. 

for  the  cavitating  hydrofoil  problea  (see  Figures  2  and  15) 
x-conponent  of  the  perturbation  velocity  with  respect  to  0., 
for  the  non-cavitating  hydrofoil  problen  (see  Figure  1) 

Free  streaa  velocity 

y-coaponent  of  the  perturbation  velocity  with  respect  to 

for  the  cavitating  hydrofoil  problen  (see  Figures  2  and  15) 

Induced  velocity  distribution  in  the  wake  of  the  hydrofoil 

y-component  of  the  perturbation  velocity  with  respect  to 

for  the  non-cavitating  hydrofoil  problen  (see  Figure  1) 

Nondinensional  v.  (x),  defined  as:  v(x)  “  v  (x) 
m 


o.D 

2 

Cavity  voluae,  non-dimens  ionalized  with  respect  to  c  ;  for 
supercavitating  hydrofoils,  includes  also  the  foil. 
Coordinate  systea  for  the  hydrofoil  (see  Figure  1) 
x-coordinate  of  the  ainiaua  pressure  point  on  the  hydrofoil 


Transforned  x  variable  defined  as  :  ■ 


partially  cavitating  and  as  z  ■  J-  *  ;  0<x<*  for  the 

I  i-x 


x  ;  0<x<l  for  the 
1-x 


supercavitating  hydrofoil  problen 
_  Angle  of  attack  of  the  free  streaa  with  respect  to  the 
hydrofoil  (see  Figure  1) 

Ideal  angle  of  attack 

Cavity  vortex  distribution  defined  as: 

Y(x)  -  Y  (x)  -Y  (x) 

c  w 

Hondiaensional  vortex  distribution  for  the  supercavitating 
hydrofoil  problea,  defined  as:  ^(x)  *  — y 

CD 

Vortex  distribution  for  the  cavitating  hydrofoil 
problea  (see  Figures  2  and  15) 
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discrete  cavity  thickness  distribution 
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F: 

q(x): 

q(x): 

qc(x) : 

qw(x): 

q#(x) : 
Q(Qc): 


r: 


R: 

R(u,*): 


s  : 
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t: 


Vapor  presure  of  the  fluid  or  pressure  inside  the  cavity 
Free  strean  pressure 
Constant  defined  in  (18.5) 

Cavity  source  distribution  defined  as  q(x)  -  qc(x)  -  qy(x) 
Hondinensional  source  distribution  for  the  super eavit at ing 

qc(*) 

hydrofoil  probles,  defined  as:  q(x)  *  — jj — 

*  ® 

Source  distribution  for  the  cavitating  hydrofoil 
problen  (see  Figures  2  and  15) 

Source  distribution  for  the  non-cavitating  hydrofoil 
problen  (see  Figure  1) 

Surface  velocity  on  the  hydrofoil  or  on  the  cavity 

Any  quantity  depending  linearly  on  the  solution  o,Y(x),q(x) 

(o,yc(x) ,qc(x))  for  the  partially  cavitating 

(supercavitating)  hydrofoil  problen 

Any  quantity  depending  linearly  on  the  solution 

Y  (x),q^(x)  for  the  non-cavitating  hydrofoil  problen 

2  1/4 

Constant  defined  as  r  ■  (1  ♦  s  )  for  the  partially 
cavitating  and  as  r  “  (1+t2^^  for  the  supercavitating 
hydrofoil  problen 
Constant  defined  in  (18.5) 

Integral  defined  in  (A. 24a) 

Constant  in  the  partially  cavitating  hydrofoil  problen  defined 
as  s  -  (*/(l-*))1/2 

1  /2 

Constant  defined  as  s  1  (&o/(l-*o))  ' 

o 

Constant  in  the  supercavitating  hydrofoil  problen  defined  as  t 
-  1 / ( £—1 ) 1 /2 
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NOMENCLATURE 


a(PL;£):  Leading  edge  radius  correction  factor  defined  in  (21.9) 
a(P.  ;*•,£  ):  Integral  defined  in  (22.21) 

Lk  O 

A:  Constant  defined  in  (5.2) 

B:  Constant  defined  in  (5.2) 

c:  Chord-length  of  the  hydrofoil 

c(x):  Camber  distribution  of  the  cavity  (see  Figure  15) 

^1 *^2*^3 *^4*^5 :  Constants  defined  in  Appendix  C 
Cpi  Drag  coefficient 

C^:  Lift  coefficient 

C^:  Moment  coefficient 

it  2 

D:  Constant  defined  as  D  - - -r.  (r  -1) 

4.r 

E:  Constant  defined  in  (6.7) 

f(n;s):  Function  defined  in  (5.9) 

f(n;s,8  ):  Function  defined  in  (22.16) 
o 

F:  Constant  defined  in  (6.8) 

h(x) :  Cavity  thickness  distribution — for  supercavitating  hydrofoils 

includes  also  the  foil  thickness  (see  Figure  15) 
h^:  Discrete  cavity  thickness  distribution 

l(n;s):  Integral  defined  in  (5.7) 

I(r';s,8o):  Integral  defined  in  (C-3) 

J(fi;t):  Integral  defined  in  (A. 23) 

K:  Number  of  uniform  intervals  to  apply  Simpson's  integration 

^  :  x  coordinate  of  the  start  of  the  cavity 

x  coordinate  of  the  end  of  the  cavity 
m;  Extent  of  the  for  NACA  a  ■  m  camber  meanline 

M:  Constant  defined  in  (7.6a) 

N;  Constant  defined  in  (7.6b).  Also  number  of  points  for  the 


The  dynamic  boundary  condition 


uc(x,0+)  =  y.u^,;  0*<x<£ 


(2.4) 


The  perturbation  velocity  flow  field  (uc,vc)  can  again  be 

generated  by  a  distribution  of  vortices  y  (x >  and  sources  q  (x)  on 

the  slit  from  x«0  to  x-1.  The  problem  can  be  again  decomposed  into 

the  effective  camber  and  effective  thickness  problems  as  shown  in 

Figure  2,  provided  that  the  cavity  thickness  distribution  is  known. 

The  boundary  conditions  (2.2),  (2.3)  and  (2.4)  in  terms  of  the 

unknown  distributions  Y  (x)  and  q  (x)  will  become: 

c  Mc 


qc(x>  !  fYc(0.d5  3nu 

vc(x,0+)  =  2  +  2 I  ~T^  U-*3T' 


(2.5) 


for:  £<x<l 


q„(x)  ,  Y  (O.d? 

v  fX,o-) - V-  +  h-  I  — 


(2.6) 


for:  0<x<l 


Y(x) 

u  (x,0+)  =  — - — 
c  ^ 


1  qn(0-dz 
1  _£ - 

0 


(2.7) 


for:  0<x< £ 


The  Kutta  condition  at  the  trailing  edge  of  the  hydrofoil  has  also  to 
be  introduced  as: 


Y  (Definite 
c 


(2.8) 


The  cavity  thickness  distribution  h(x)  can  be  determined  from  the 


kinematic  boundary  condition  on  the  cavity  surface: 


But 


vc(x,0+)  = 


n  (x)  and 
u 


qc(x)  ,  1  ( 1  „  3 

2  2k  '  f  C-x  '  U-'3^ 

for:  0<x<  £ 

n^(x)  can  be  expressed  as: 


(Vh) 


(2.9) 


nu(x)  *  n(x)  +  — ^  (2.10) 

n^(x)  =  n(x)  -  (2.11) 


Subtracting  (2.6)  from  (2.9)  and  using  (2.10),  (2.11)  and  (1.6),  we 
get: 


qc(x)  -  qw(x)  =  0<x<£ 


(2.12) 


Substituting  (2.10)  and  (2.11)  in  (2.5)  and  (2.6)  and  having  in  mind 
(1.1)  and  (1.6),  we  get: 


qc(x)_qw(x)  1  r1  Yc(C)_Yw(C) 
?  +  2tt  '  T  - 


0 


C-x 


•  dC  -  0;  £<x<l  (2.13) 


qc(x)-qw(*)  ,  1  YC(S)-YW(» 

2  2.  *  $ 


C-x 


dC  =  0;  0<x< 1  (2.14) 


Combining  (2.13)  and  (2.14)  we  finally  get: 


q  (x)  =  q  (x) ;  £<x<l 


(2.15) 


Vx)-qy(x)  1  YC(€)-YW(S) 

9  J  F— v 


-  •  <U  -  0;  0<x<l  (2.16) 


We  define  the  cavity  source  q(x)  and  cavity  vortex  y(x) 
distributions  as: 


q(x)  =  q  (x)  -  qw(x) 


y(x)  -  y  (x)  -  y  (x) 
c  w 


Then  the  equations  (2.12),  (2.15)  and  (2.16)  become: 


(2.17) 

(2.18) 


q(x)  =  U  .  — ;  0<x<£ 
0°  3x 


(2.19) 


q (x)  =  0;  £<x<l 


(2.20) 


_  +  1  .  f  xL Lhik  .  0;  0<x<l  (2-21) 


The  dynamic  boundary  condition  (2.7)  can  also  be  expressed  in  terms 
of  q(x) ,  y(x)  as: 

xSsi  1_  v  ..V!i  +  JL 

9  *  7  f-v  2  '  00  2  2lT  '  Z  E~X 


2  2tt  '  C-X  2  00  2  2tt  ^  £-x 


(2.22) 


0<x<£ 


The  right  hand  side  of  (2.22)  is  the  negative  of  the  perturbation 
velocity  uv(x,0+)  on  the  upper  side  of  the  hydrofoil  in 
non-cavitating  flov,  as  it  can  be  seen  from  (1.1)  and  (1.7). 


The  Kutta  condition  will  become: 


Y(l)-finite  (2.23) 

In  the  preceding  analysis  nothing  has  been  said  about  the  trailing 

edge  of  the  cavity.  In  the  nonlinear  theory  using  Riabouchinsky  or 

reentrant  jet  models  it  has  been  proven  in  [13]  that  the  thickness  at 

2 

the  end  of  the  cavity  is  of  order  of  c  .  Therefore  in  linear 
theory  we  require  that: 

h(£)  =  0  (2.24) 

or  by  integrating  (2.19)  we  get  the  so-called  cavity  closure 

condition:  £ 

/  q(C).dC  =  0 

0  (2.25) 

The  final  system  of  singular  integral  equations  is: 


q(x) 

2 


+  -k-  f 
0 


Y(p,d^  =  0;  0<x<l 
(-x 


(2.21) 


y(x)  1 

2  2tt 


£ 


i 

0 


qCQ.dg  o 
C-x  r  2  ' 


U  =  -  u  (x, 0+) 

oo  W 


for:  0<x<£ 
y(1)  ■  finite 


£ 

/  q(0.dC  =  0 
0 


(2.26) 

(2.23) 

(2.25) 


The  unknowns  are  a,  q(x)  and  y(x).  The  solution  corresponding  to 

a  certain  cavity  length  £  depends  only  on  the  distribution  of 

u  (x.0+)  for  0<x< £  and  is  unique,  as  will  be  proven  in  Section  4. 

* 
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The  non-cavitating  hydrofoil  problem  can  be  decomposed,  in  linear 
theory,  into  the  camber,  the  thicknesa  and  the  angle  of  attack  (or  flat 
plate)  problems.  If  is  any  quantity  depending  linearly  on  the 
distributions  q^(x)  and  Y^(x),  then: 

.  qw-C  +  qwt  + C 

vhere  the  superscripts  c,T,a  correspond  to  the  camber,  the  thickness 
and  the  angle  of  attack  problems  respectively.  The  total  camber  n(x) 
has  been  decomposed  as: 

ri(x)  =  n(x)  +  a.(l-x)  (3.2) 

where  n(x)  is  the  hydrofoil  camber  and  a  the  angle  of  attack  as 
defined  in  Figure  1. 

If  we  put  Q  «uw(x,0+)  in  (3.1)  then: 

(3  3 ) 

u  (x.CH-)  =  u  C(x,0+)  +  u  T(x,0+)  +  u“(x,0+) 

W  w  w 

• 

Since  the  right  hand  side  of  equations  (2.21),  (2.26),  (2.23),  (2.25) 
can  be  decomposed  as  in  (3.3),  the  same  is  true  for  the  solution  a, 
q(x) ,  Y (x) : 

o  =  aC  +  aT  +  aa  (3.4) 

q  =  qC  +  qT  +  q°  (3.5) 

c  .  x  ,  a 

Y  =  Y  +  Y  +  Y  (3.6) 

Furthermore,  for  any  quantity  Q  depending  linearly  on  the  solution  it 
will  be  true: 

Q  -  QC  +  QT  +  Q°  (3.7) 

where  Q  can  be,  for  example: 

-the  cavity  thickness  distribution  h(x) 


-the  cavity  volume  V 
-the  lift  coefficient, 

-the  moment  coefficient, 

-the  square  root  of  drag  coefficient. 

Therefore,  the  partially  cavitating  hydrofoil  problem  can  be 
decomposed  into  the  camber,  the  thickness  and  the  flat  plate  problems, 
all  corresponding  to  the  same  cavity  length  £,  as  it  is  depicted  in 
Figure  3.  Some  of  these  elementary  problems  possibly  have  no  physical 
meaning,  i.e.  they  result  in  negative  or  convex  cavities,  but  if 
combined  with  the  appropriate  angle  of  attack  they  give  meaningful 
solutions  to  our  problems  as  illustrated  in  Figures  10  and  11. 

As  an  application,  the  following  special  cases  can  be  examined: 

a)  Change  in  angle  of  attack: 

If  we  change  the  angle  of  attack  of  a  partially  cavitating 
hydrofoil  by  Aa,  then  the  solution  for  the  same  cavity  length 
will  change  by  an  amount  which  is  the  solution  to  a  flat  plate  at  angle 
of  attack  Aa  and  the  same  cavity  length. 

b)  Change  of  the  camber/chord  ratio: 

For  a  given  camber  meanline  shape  partially  cavitating  at  zero 
angle  of  attack,  if  we  change  the  camber /chord  ratio  by  a  factor,  then 
the  solution  for  the  same  cavity  length  changes  by  the  same  factor. 

A  similar  relationship  holds  for  affine  thickness  forms. 

Therefore  for  a  given  camber  or  thickness  series  partially 
cavitating  at  a  cavity  length  £,  we  need  to  solve  the  problems  for  one 
value  of  the  similarity  parameter  and  one  angle  of  attack. 


Partially  cavitating 
hydrofoil  problem 


| 


Figure  3.  Decomposition  of  the  partially  cavitating  hydrofoil  problem 
into  the  camber,  thickness,  and  flat  plate  problems. 


We  have  to  solve  the  linear  system  of  singular  integral  equations 
(2.21),  (2.26)  in  combination  with  the  restrictions  (2.23),  (2.25). 

We  first  invert  the  equation  (2.21)  by  using  the  formulas  (1.1)  and 
(1.3): 


Y(x) 


dC 


(4.1) 


The  solution  (4.1)  already  satisfies  the  Kutta  condition  (2.23)  and  is 
unique  as  long  as  q(x)  is.  Substituting  in  (2.26)  y(x)  from  (4.1) 
we  get  the  following  integral  equation  for  q(x): 


h  { <4-2> 

We  now  introduce  the  following  transformation  as  suggested  by  Hanaoka 
in  [111: 


"  \  1-x  ’  n  "  ^Jl-C  ’  S  "  {l-i 


0<x,  £<Sl  and  0<z,  n<s 


(4.3) 


and  equation  (4.2)  becomes: 


o.U 


i_  i  n^alnl  -dn_  =  I  _ z .  ( — •  _  u  +>  (4  4) 

2tt  ’  *  2  '  (z-n)  2  2  V  2  Uw  ' 

0  (1+n  )  v  n;  ,-L- 


1+z 


which  is  in  a  form  similar  to  equation  (1.1). 


Inversion  of  (4.4)  in  combination  with  the  requirement  that  the 


cavity  source  is  square  root  singular  at  the  trailing  edge  of  the 
cavity  [43],  renders  the  solution. 


q(z) 


_1 

IT 


(1+z2)  rs  Vn(s-n) 

\/z.  (s-z)  0  1+n2 


1 

z-n 


dn  (4.5) 


which  is  unique  as  long  as  o  is. 

Equations  (4.5)  and  (4.1)  give  us  the  formal  unique  solution  to 
the  general  partially  cavitating  hydrofoil  problem  as  soon  as  we 
determine  the  corresponding  cavitation  number  a  • 


5.  Cavitation  number  n 

By  using  (A. 19)  equation  (4.5)  becomes: 

■  ^ toll  .  u  +(n)  .  - 

^{s-z)  l  1+f)2  w  z-n 


q<z)  =  7 


1 


dn  +  o.U 


A-B.  z 


00  *  ^z. (s-z) 


(5.1) 


where : 


A  •  JT?  •  B  =  T^f  ’  f**'’  r2  =  7PT  (5-2> 

To  determine  o  we  have  to  apply  the  closure  condition  (2.25) 
which  in  terms  of  the  transformation  (4.3)  becomes: 

JSq(z)  .  -  0  (5.3) 

0  (1+Z  ) 

Substituting  (5.1)  in  (5.3)  and  using  (A. 11)  and  (A. 12)  we  get: 


2 _  rs  z  _ 1 _  ,  /  Vn(s-n)  +  _J_ 

°  '  U.-D  '  0  *'/2(s'2)  ‘  l+z2  '  0  l+„2  "  ‘  Z'n 


(5.4) 


where : 


D  =  •  (rZ-l) 

4.r4 


(5.5) 


By  using  the  Poincare-Bertrand's  formula  (A. 28)  >we  change  the  order  of 
integration  in  the  double  integral  of  (5.4)  and  we  get: 


-U S  i  ^  •  C  •«"")•  0" 


(5.6) 


where : 


Kn;s)  =  j  •  dz 

0  Vs  2  i+z2  2  n 


(5.7) 


The  singular  integral  I(  n;s)  is  evaluated  analytically  in  Appendix  A. 
Thus  the  final  expression  for  a  becomes: 


a  2.u  (n) 

a  =  /  - -  .  f(n;s)  .  dn 


(5.8) 


where : 


.JB5I  -|L. 

d+n  V  r  -1  (5.9) 


with 


r2  =\h+s2 


(5.10) 


Some  plots  of  the  function  f(  n  ;s)  for  different  values  of  l  can  be 
seen  in  Figure  6. 

In  other  words,  a  and  q(z)  have  been  expressed  in 


• »  .  \  %  %  .  %  ,  *  .  *  %  .  •  A  N.  i*,  \ 


(5.1)  and  (5.8)  in  terms  of  integrals  of  known  functions  over  the 
cavity  length.  We  observe  again  the  dependence  of  the  solution  only  on 
the  distribution  of  the  upper  surface  perturbation  velocity  u^+  in 
the  non-cavitating  flow.  Similar  expressions  can  also  be  derived 
for  y(x),  CD  (see  also  [16].) 

In  this  work  we  focus  more  on  the  cavity  shape  than  on  the 
hydrodynamic  characteristics  of  the  cavitating  hydrofoil.  Thus  we 
proceed  by  giving  general  expressions  for  the  cavity  volume  V  and  the 
thickness  distribution  h(x). 


6.  Cavity  Volume  V 

The  cavity  thickness  distribution  comes  from  the  integration  of 
(2.19): 


(6.1) 


where  we  have  made  use  of  the  boundary  condition  at  the  leading  edge  of 
the  cavity: 

h(0)-0  (6.2) 


The  cavity  volume  V  will  be  given  as: 


V  =  /  h(x).dx  =  x.hl  -  /  x.-^.dx  =  ”  /  x'^u  ,c*x  (6*3) 


0  0 


where  we  have  made  use  of  (6.2),  (2.24)  and  (2.19). 


In  terms  of  the  transformation  (4.3)  equation  (6.1)  becomes: 


which  by  using  (4.5),  (A.28)  and  (A.21)  finally  renders: 


where: 


(6.5) 

(6.6) 

(6.7) 

(6.8) 
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7.  Analytical  solutions  for  some  hydrofoils 

As  a  first  step  to  our  solution  we  have  to  solve  the  non-cavitating 
hydrofoil  problem  in  order  to  find  u+  ,  given  as: 


u„+<x) 


Y.  (x) 


1  q  (O-dC 


=  u  (x,0+)  =  — « - •  7  r 


(7.1) 


For  the  most  camber  meanlines  used  in  practice,  the  vortex  distribution 
Y^(x)  is  known,  since  they  are  designed  to  develop  that  load. 

The  thickness  source  distribution  q^(x)  is  related  to  the 
thickness  distribution  t(x)  as  in  (1.6).  The  integral  in  (7.1) 
corresponds  to  the  induced  perturbation  velocity  u^  (x,0)  due  to 
the  foil  thickness,  and  for  some  special  cases  can  be  evaluated 
analytically. 

Analytical  results  to  the  partially  cavitating  hydrofoil  problem 
can  be  obtained  as  long  as  the  integral  in  the  expressions  (5.1),  (5.8) 
and  (6.5)  can  be  evaluated  analytically.  The  following  cases  are 
considered  as  illustrated  in  Figure  4. 

a.  Flat  plate 

For  the  non-cavitating  flow  we  have  by  using  (1.3)  with  =  -a: 


and  by  using  (7.1)  and  (4.3) 


Yw(x)  =  2-U^.a  0<x<l 


(7.2) 


Uoo*a 

u  +(n)  =  — ;  0<n<_s 

w  n 


(7.3) 


Substituting  (7.3)  in  (5.8)  and  (5.1)  and  by  using  the  formulas  from 
Appendix  A,  we  get  (see  also  [10]  and  [3l]): 
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(d)  Planoconvex  hydrofoil 


Figure  4.  Special  cases  of  partially  cavitating  hydrofoils 


Table  5 


CAVITATION  NUMBER  a  FOR 
PARABOLIC  CAMBER  OF  X  =  0. 4  AT  a  =  0° 


Numerical 

l 

Analytical 

K  =  20 

K  =  10 

K  =  5 

0.  1 

0.4929822 

0.4929822 

0.4929822 

0.4929822 

0.3 

0.8047841 

0.8047841 

0.8047841 

0.8047825 

0.5 

0.9656854 

0.9656854 

0.9656854 

0.9656083 

0.7 

1.0359341 

1.0359341 

1.0359341 

1.0362429 

0.9 

0.9989466 

0.9989466 

0.9989687 

1.0054026 

Table  6 


CAVITY  VOLUME  V  FOR 
PARABOLIC  CAMBER  OF  A  =  0.4  AT  a  =  0° 


Numerical 

l 

Analytical 

K  =  20 

K  =  10 

- - 

K  =  5 

0.  1 

-0.0003569 

-0.0003569 

-0.0003569 

-0.0003569 

0.3 

-0.0050591 

-0.0050591 

-0.0050591 

-0.0050595 

0.5 

-0.0162661 

-0.0162661 

-0.0162661 

-0.0162791 

0.7 

-0.0332860 

-0.0332860 

-0.0332860 

-0.0331346 

0.9 

-0.0550233 

-0.0550233 

-0.0550161 

-0.0554519 

Table  4 


FOR 

CAVITY  THICKNESS  DISTRIBUTION  h(x) 

FLAT  PLATE  AT  a  =  1  RAD  AND  £  =  0.5 

X 

Numerical 

N  «  10 

N  *  20 

N  =  40 

0. 

0. 

0. 

0. 

0.0006 

0.0188053 

0.0188283 

0.0188298 

0.0090 

0. 1370588 

0. 1371007 

0. 1371033 

0.0407 

0.3877945 

0.3878217 

0.3878234 

0. 1066 

0.6975954 

0.6975761 

0.6975750 

0.2 

0.9307760 

0.9307408 

0.9307387 

0.3 

0.9944224 

0.9944067 

0.9944058 

0.3866 

0.8845249 

0.8845262 

0.8845263 

0.45 

0.6487553 

0.6487611 

0.6487616 

0.4876 

0.3402218 

0.3402258 

0.3402261 

0.5 

-0.0000009 

-0.0000002 

-0.0000001 
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CAVITY  SOURCE  DISTRIBUTION  q(x) 

FOR  FLAT  PLATE  AT  a  -  1  RAD  AND  t  ■*  0.5 
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Table  1 


a/°  FOR  FLAT  PLATE 

l 

Analytical 

Numerical 

K  =  20 

K  -  10 

K  =  5 

0.1 

0.0395011 

0.0395011 

0.0395011 

0.0395011 

0.3 

0.0679238 

0.0679238 

0.0679238 

0.0679238 

0.5 

0.0857864 

0.0857864 

0.0857864 

0.0857862 

0.7 

0.0956519 

0.0956519 

0.0956519 

0.0956511 

0.9 

0.0865823 

0.0865823 

0.0865823 

0.0866735 

Table  2 


CAVITY  VOLUME  V  FOR  FLAT  PLATE  AT  a  =  1  RAD 

l 

Analytical 

Numerical 

K  =  20 

K  =  10 

K  =  5 

0.1 

0.0313815 

0.0313815 

0.0313815 

0.0313815 

0.3 

0.1676164 

0.1676163 

0.1676164 

0. 1676189 

0.5 

0.3758550 

0.3758550 

0.3758549 

0.3758483 

0.7 

0.6699281 

0.6699281 

0.6699281 

0.6690742 

0.9 

1.2007285 

1.2007285 

1.2007252 

1.2203253 

now  observe  that  for  the  ine  flow  conditions  the  linesr  theory 
overnredicts  the  extent  end  the  volume  of  the  cavity  for  «  NACA  16-006 


hydrofoil. 


my  t«  r.»  f.'vjvsr.  r.  r;v. 


l1*  :  •  n  •  i-  «-  »■■ 


increase  the  thickness  and  keep  the  same  flow  conditions,  the  cavity 
becomes  longer  and  has  more  volume.  This  is  the  reverse  from  what  we 
concluded  for  biconvex  hydrofoils  and  is  due  to  the  change  of  camber 
for  a  planoconvex  hydrofoil  as  we  increase  its  thickness.  It  also  is 
the  reverse  from  what  we  observe  in  the  experiments  [26,  29,  and  37]. 

Some  cavity  plots  for  a  biconvex  and  a  planoconvex  hydrofoil  and 
for  each  of  the  elementary  problems  as  described  in  Section  3  can  be 
seen  in  Figures  10  and  11.  We  observe  that  the  thickness  and  camber 
problems  are  physically  unaccepted,  since  they  result  in  negative 
cavities  or  even  negative  cavitation  numbers.  Nevertheless,  those 
problems,  combined  with  some  angle  of  attack,  give  meaningful  solutions 
to  our  problem. 


10.  Comparison  of  the  linear  to  the  nonlinear  theory 

The  linear  theory  described  in  the  previous  sections  is  compared 
with  the  nonlinear  theory  using  the  modified  Riabouchinsky  cavity 
termination  model. 

First,  for  a  flat  plate,  our  results  and  the  nonlinear  results 
obtained  by  using  the  analysis  described  in  [33]  are  shown  in  Figures 
12  and  13  for  a  m  2°  and  a  *  5°  respectively.  We  can  anticipate  that 
for  the  same  flow  conditions  the  linear  theory  underpredicts  the 
extent  and  the  volume  of  the  cavity  for  a  partially  cavitating  flat 
plate.  We  also  observe  that  the  discrepancies  between  linear  and 
nonlinear  theory  become  smaller  as  we  decrease  the  angle  of  attack. 

The  cavity  plots  obtained  from  linear  and  nonlinear  theory  [38] 
for  a  NACA  16-006  thickness  form  at  a  ■  4°  are  shown  in  Figure  14.  We 
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The  numerical  integrations  are  performed  for  a  flat  plate  and  a 
parabolic  camber,  and  the  results  are  compared  vith  the  analytical  ones 
described  in  Section  7.  The  agreement  is  excellent,  even  for  a  small 
number  of  Simpson's  intervals  X  and  for  the  whole  range  of  cavity 
lengths  as  shown  in  Tables  1,  2,  3,  5,  and  6.  The  cavity  thickness 
distribution  for  a  flat  plate  eavitating  at  lm  0.5  is  computed  by 
using  the  algorithm  described  by  the  expressions  (8.5),  (8.6)  and 
(8.7).  The  results  are  shown  in  Table  4  for  a  different  number  of 
points  H,  where  K  ■  50  has  been  used  for  the  evaluation  of  the  cavity 
source  distribution  q(x) .  The  results  seem  to  be  quite  insensitive  to 
H,  since  a  four-decimal-places  accuracy  is  obtained  rather  quickly. 

9.  The  biconvex  and  the  planoconvex  hydrofoils 

The  numerical  computations  as  described  in  the  previous  section 
are  carried  out  for  a  biconvex  and  a  planoconvex  hydrofoil.  In  Figures 
7,  8,  and  9,  graphs  of  the  results  are  shown  for  4Z  and  8Z  maximum 
thickness/chord  ratios.  From  those  graphs  the  following  conclusions 
can  be  driven: 

a)  A  biconvex  hydrofoil  seems  to  simulate  a  flat  plate  better 
than  a  planoconvex  hydrofoil  vith  the  same  thickness,  even  though  the 
last  two  have  the  same  lower  surfaces. 

b)  For  biconvex  hydrofoils  and  short  cavities,  if  we  increase 
the  thickness  and  keep  the  same  flow  conditions,  the  cavity  becomes 
shorter  and  has  less  volume. 

c)  For  planoconvex  hydrofoils  and  short  cavities,  if  we 


is  singular  and  needs  more  careful  numerical  treatment.  Instead,  by 
using  (A.  19)  ve  put  (5.1)  in  the  form: 


q(z) 


_  +.  .  + 

1+z2  fSVn(s-ri)  Uw  Z  uw 

11  *  izCs’z)  '  o  l+n2  *  z'n' 

+  t°-u0O-2'uw  (z)1  •  yz(s-z) 


dn  + 


(8.3) 


u  (z)  -  u  (n)  du  (n) 
The  integral  in  (8.3)  is  not  singular  since  — - - - ^  — 2— - — 

z-n  d(n) 

as  rr*-z  and  therefore  it  is  evaluated  by  applying  Simpson's 
rule  in  6  . 

To  find  the  cavity  thickness  distribution  ve  perform  the 
integration  in  (6.1)  numerically  after  making  the  transformation: 

5  *  J,.sin2(-|) 

0<£<£  0<_4>^_n  (8.4) 

We  then  compute  the  discretized  thickness  distribution  by  using  the 
algorithm: 


i+1 


i+1 


hi  +  2fiT 


q(C)  - s±nd> . d<p ;  i  *  1,N 


with: 


where : 


«  0  at  <t>i  =  0 


i+1 


1  =  1,N 


(8.5) 

(8.6) 

(8.7) 


The  integral  in  (8.5)  is  computed  numerically  by  applying  Simpson's 
rule  with  one  interval  (i.e.  three  points).  A  check  to  the  numerical 
integrations  should  be  the  closure  condition: 


(7.22) 


uw+(x)  *  e.U— ;  0<x<l 
Therefore  the  partially  cavitating  problem  ia  identical  to  the 
non-cavitating  one,  once  ve  put: 

o  -  2. c  (7.23) 

Again  that  reaulta  to  aero  cavity  thickneaal 

8.  Wuaerical  method 

For  general  shape  hydrofoils,  u*(x)  can  be  evaluated  either 
analytically  or  numerically.  Then  the  integrals  involved  in  the 
formulas  (5.8),  (5.1)  and  (6.5)  will  be  evaluated  numerically. 

Ve  first  define  the  transformation  variable  8  as  follows: 


4  2/0N 

n  *  s.sin  Or) 


0<n<s  -*-*•  0<8<tt 


(8.1) 


Then  ve  perform  the  numerical  integrations  in  terms  of  6,  by  applying 
Simpson's  rule  with  K  uniform  intervals  (i.e.,  2K+1  points). 

For  example  equation  (5.8)  becomes: 

+  o 

sin  6 _ _ 

(8.2) 


2  2  it 

a  -  —  .  V1"  •  /  <U^>  * 

r  -1  0  - 


4L 


(l+s2.sin4(|))2  ‘ 


[^+1  -  s.sin2(|)  .  ^r2-l]  ,d6 

The  integrand  in  (8.2)  behaves  smoothly  at  8«0  and  8 ■it  and 

therefore  ve  expect  more  accurate  results  by  applying  Simpson's  rule  in 

(8.2)  than  in  (5.8)  in  which  the  integrand  behaves  like  square  root  of 

the  distance  from  the  ends.  Furthermore  for  the  case  u+  behaves  like 

v 

1 h  at  n-0  (as  in  (7.3))  then  by  using  (8.2)  ve  avoid  the  square  root 

singularity  at  n-0  in  the  integrand  of  (5.8).  The  integral  in  (5.1) 
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Y  (x)  -  Vi  0<x<« 
w  1 


as: 


and: 


(7.16) 


u  +(x)  ■  A.U  ;  0<x<m 

y 

u  +(x)  -  A.U  .  (4=5)  :  «<x<l 
w  2  •  1ID 


(7.17) 


For  £<m  the  non-cavitat ing  boundary  value  problem,  aa  ve  can  see  in 
Figure  5a,  already  satisfies  the  partially  eavitating  problem  stated  by 
the  equations  (2.21),  (2.26),  (2.23)  and  (2.25),  if  we  put: 


c  ’ Yi . 

(7.18) 

The  respective  unique  solution  will  be: 

q(x)  -  0 

(7.19) 

Y(x)  •  0 

(7.20) 

In  other  words  there  is  no  cavity  for  £<*! 

The  above  results  can  also  be  obtained  by  substituting  (7.17)  in 
(5.8),  (5.1)  and  (4.1).  The  equation  (5.8)  in  combination  with  (7.17) 
and  (7.18)  requires  the  function  f(n,a)  to  satisfy  the  restriction: 

s 

/  f(n;s).dn  “  1 

0  (7.21) 
which  also  can  be  proved  from  direct  integration.  In  other  words,  the 
function  f(n;s)  becomes  a  "delta"  function  as  £-*0,  as  can  also  be 
seen  in  Figure  6. 

f.  Ellipse  at  a“0  degrees 

Consider  an  ellipse  of  axes  (l,c)  subject  to  a  uniform  flow 
V  ,  as  shown  at  figure  5b.  Then  according  to  linear  theory  [39]: 


a  =  2X.\/r.(l+  yi-«-  ) 

V  -  -  ^.x.d3/2.  (1-  ) .  ( 3—  2 £. ) 

c.  Biconvex  hydrofoil 
This  is  a  symmetric  hydrofoil  confined  between  two  circular  arcs. 

The  thickness  distribution  to  first  order  is  given  as: 

t(x)  =  4.u.x. (1-x) 

By  inserting  (7.13)  in  (1.6)  and  (7.1)  we  get: 

.  2.yU 

u  (x)  =  — -  .  [2  +  (l-2x)  .  *!>(-£-)];  0<x<l 

W  7T  i  X 

or  in  terms  of  the  transformed  variable  n 

4.p.U  ,2 

u  +  (n)  - - 1  .  (1  +  i-iL.  .  Jin(n));  0<n<s  (7.15) 

w  *  1+r,2 

The  integrals  involved  in  the  expressions  (5.8)  and  (6.5)  are  computed 
numerically  as  described  in  Section  8. 

d.  Planoconvex  hydrofoil 

This  is  the  hydrofoil  with  the  upper  side  a  circular  arc  and  the 
lower  side  flat.  It  can  be  considered  as  a  combination  of  a  parabolic 
camber  and  a  biconvex  hydrofoil  with  w  ■  A  as  illustrated  in 
Figure  4.  Therefore  the  results  will  be  the  superposition  of  the  two 
previous  cases,  according  to  Section  3. 

e.  BACA  a-m  camber  meanlines 

The  vortex  distribution  at  ideal  angle  a«a.  will  be  given 


(7.13) 


(7.14) 


W .11/ 

(7.12) 


q(z)  =  a'u«  *  WlTT-h"'  0<z<s 


(7.5) 


M  =  77  [~-  (r2+l)^  +  \(r2-l)h] 

N  =77 


(7.6a) 


(7.6b) 


By  using  (6.5)  and  (7.3)  we  also  recover  after  lengthy  algebra  the 
result  by  Geurst  for  the  cavity  volume  V  [10]: 


v  =  .  *(1+4.  Vl^) 

16 


(7.7) 


b.  Parabolic  camber 


In  that  case  the  camber  distribution  is  given  as: 


i(x)  =  2.A.x.(l-x);  0<x<l 


(7.8) 


Substituting  (7.8)  in  (1.3)  we  get: 


Y  (x)  =  8.U  .A.  yx(l-x);  0<x<l 
w  00  » 


(7.9) 


and  thus: 


u  +(n)  =  4.U  .A. — 0<n<s 


(7.10) 


Then  the  formulas  (5.8)  and  (6.5)  will  recover  the  results  by  Geurst 

for  a«0  degrees  [12]: 
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Biconvex  8%  thick 


ot/o  vs.  %  for  biconvex  hydrofoils,  planoconvex  hydrofoils, 
and  a  flat  plate,  all  at  a  -  4°  -  Linear  theory 


Biconvex  8%  thick 


Figure  8.  V  vs.  £  for  biconvex  hydrofoils,  planoconvex  hydrofoils, 
and  a  flat  plate,  all  at  a  =  4  -  Linear  theory 
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Biconvex  8%  thick 


»  9.  V  vs.  a/o  for  biconvex  hydrofoils,  planoconvex  hydrofoils, 
and  a  flat  plate,  all  at  a  =  4  -  Linear  theory 


1.5 

1. 

• 

‘  0°,  a  =  0.0966 

(b)  Biconvex  hydrofoil  at  a  *  0  ,  o  *  -0.0013 


(c)  Flat  plate  at  a 


0.8138 


(d)  Planoconvex  hydrofoil  at  a  ■  4  ,  a  *  0.9091 

Figure  11.  The  elementary  problems  for  partially  cavitating  4%  thick 
planoconvex  hydrofoil  at  Z  *  0.5  (a)  +  (b)  +  (c)  *  (d) 


V-:  v'.-.'  -.'.-.' --V-V .V«V*VAV>V-V -V*' -V-’v'-V.;  .  *.‘A 


• 

0.4 

0.7 

(a)  “nonlinear  *  °'659’  “linear  ’  °'623 

nonlinear  0.424,  °iinear 


(c)  "nonlinear  *  °‘384’  "linear  =  °'365 

Figure  12.  Partially  cavitating  flat  plate  at  a  *  2° 

■  —  -  nonlinear  theory  [33  ] 

-  linear  theory  for  the  sr.ne  cavity  length  as  nonlinear  theory 

—  —  —  —  linear  theory  for  the  s?ne  cavitation  number  as  nonlinear  theory 


V.*’ 


0.1 


(c)  a 


nonlinear 


1.028,  o 


linear 


0.912 


Figure  13.  Partially  cavitating  flat  plate  at  a  *  5 
i  —  — —  nonlinear  theory  [  33  ] 

-  linear  theory  for  the  same  cavity  length  as  nonlinear  theory 

-  linear  theory  for  the  same  cavitation  number  as  nonlinear  theory 


The  condition  at  infinity: 


u  -*0  and  v  -*-0  at  infinity 
c  c 


(11.6) 


Thn  boundary  conditions  in  tarns  of  the  unknown  distributions 


r  (x)  and  q  (x)  will  bacons: 
c  c 


Yc(x)  1  f  0 

2 - f  ~~"x  -  2  *  U-  "  0  :  0  <  x  "  1 

0 

with  y  (x)  -  0  for  1  <  x  <  l 


(11.7) 


«e,*)  .  1  Tc<*>'d« 

2  2.it  *  5  -  x 

o 


-  3x 


0  <  x  <  1 


Yc(l)  -  0 


/  qc(x)  •  dx  -  0 
o 


(11.8) 


(11.9) 

(11.10) 


The  boundary  condition  (11.6)  is  autonatically  satisfied  by  the  use  of 

vortices  and  sources  for  the  representation  of  the  velocity  field.  For 

liven  1  and  n£(x),  the  equations  (11.7)  -  (11.10)  give  us  the  unique 

solution  o  ,y  (x)  and  q  (x) . 

c  c 

The  cavity  thickness  h(x)  and  the  cavity  caaber  c(x)  will  then  be 


given  by  the  following  relationships: 

•  Hr  •  ’<>>  i  o  <  *  <  t 


(11.11) 


c(x)  -  n£(x)  + 


;  0  <  x  <  1 


(11.12a) 


U-  '  H  •  ’in00  *  1  <  *  *  ‘ 


(11.12b) 


,<*>  *  -  fc  •  / 


i  !  V5)-d« 


x  -  C 


;  1  <  x  <  l 


(11.13) 


where  v^fl(x)  the  induced  velocity  in  the  wake  by  the  vorticity 

distribution  7  (x) . 

c 


12.  Decomposition  in  the  supercavitating  hydrofoil  problem 


If  the  hydrofoil  has  a  camber  meanline  n(x),  a  thickness 
distribution  t(x)  and  is  at  an  angle  of  attack  a ,  as  shown  in  Figure 
16,  then  its  lower  surface  can  be  expressed  as: 


and 


n0  (x)  =  n(x)  -  +  a  •  (l-x) 


(12.1) 


9n£  =  9n  _  1  .  3t 
3x  9x  2  9x 


(12.2) 


By  decomposing  the  RHS  of  the  system  of  equations  (11.7)  -  (11.10) 
as  in  (12.2),  we  can  easily  see  that  as  long  as  we  keep  £  fixed,  the 
supercavitating  general  shape  hydrofoil  problem  can  be  considered  as 
the  superposition  of  the  following  three  elementary  problems  (also 
illustrated  in  Figure  16) : 

The  camber  problem  with  n^(x)  =  0 (x) 

_  ,  ^ (x) 

The  thickness  problem  with  n^x)  - 2 — 

The  flat  plate  problem  with  n^(x)  =  <**(l-x) 

The  solution  o,  Yc(x),q»c(x)  will  be  the  linear  superposition  of  the 
solutions  to  those  problems: 


a  =  aC  +  aT  +  aa 


(12.3) 


Yc  (x) 


Yc(x)  +  Yc(x)  +  Yc(x> 


(12.4) 


qc  (x) 


*ip(x)  +  q*(x)  +  q“(x) 


(12.5) 


Furthermore,  for  each  quantity  Qc  which  depends  linearly  on  the 
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y 


ire  16.  Decomposition  of  the  S 


6‘ 


ilution,  we  will  have: 


Q  =  QC  +  QT  +  Q° 
c  c  x  c 


(12.6) 


r  example,  Q  can  be:  a,  y  (x),  q  (x),  h(x),  c(x),  V,  C  , 

C  C  C  L 

[•  Vv  Some  particular  cases  can  then  be  examined: 

Change  in  the  angle  of  attack.  If  we  change  the  angle  of 
tack  of  a  supercavitating  hydrofoil  by  Aa ,  then  the  solution  for  the 
me  cavity  length  will  change  by  an  amount  which  is  the  solution  to  a 
percavitating  flat  plate  at  angle  of  attack  Aa  and  the  same  cavity 
ngth,  as  also  illustrated  in  Figure  17a. 

Change  of  the  camber/chord  ratio.  For  a  given  camber  meanline 
ape  supercavitating  at  zero  angle  of  attack,  if  we  change  the 
mber/chord  ratio  by  a  factor,  then  the  solution  for  the  same  cavity 
ngth  changes  by  the  same  factor,  as  also  illustrated  in  Figure  17b. 
similar  relationship  is  also  true  for  affine  thickness  forms. 

For  a  given  camber  or  thickness  series,  therefore,  we  need  to 
lve  the  problem  for  one  value  of  the  similarity  parameter  and  one 
gle  of  attack. 

The  decomposition  described  above  can  also  be  done  when  the  cavity 
arts  after  the  leading  edge  of  the  foil.  In  that  case,  both  the 
ading  edge  and  the  trailing  edge  of  the  cavity  have  to  remain  fixed 
r  each  elementary  problem. 


.  Solution  to  the  general  supercavitating  hydrofoil  problem 
It  is  convenient  to  introduce  the  following  nondimens ional 
antit ies : 
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19.  o  vs.  i  for  some  Supercavitating  Hydrofoils 


Table  11 


ELLIPTIC  FOIL  WITH  e  =  1  AT  Z  =  2. 
CAVITY  THICKNESS  DISTRIBUTION  h(x) 


Numerical 

X 

Analytical 

N  -  10 

N  =  20 

0.0000 

0.000000 

0.000000 

0.000000 

0.0489 

0.618034 

0.618036 

0.618034 

0.1910 

1.175571 

1.175574 

1.175571 

0.4122 

1.618034 

1.618039 

1.618034 

0.6910  • 

1.902113 

1.902119 

1.902113 

1.0000 

2.000000 

2.000007 

2.000000 

1 . 3090 

1.902113 

1.902119 

1.902113 

1.5878 

1,618034 

1.618039 

1.618034 

1 . 8090 

1.175571 

1.175574 

1.175571 

1.9510 

0.618034 

0.618036 

0.618034 

2.0000 

0.000000 

-0.000000 

-0.000000 

Table  9 


3.14159265 


3.14159265 


3.14159217 


3.13764913 


Table  10 


■ELLIPTIC  FOIL  WITH  e  *  1  AT  £  «  2 
CAVITY  SOURCE  DISTRIBUTION  q(x) 


Analytical 


2.06474 


1.33333 


0.75000 


0.57735 


0.31449 


0.10050 


0. 


Numerical 

K  =  20 

K  =  10 

K  =  3 

2.06474 

2.06474 

2.06206 

1.33333 

1.33333 

1.33270 

0.75000 

0.75000 

0.74987 

0.57735 

0.57735 

0.57727 

0.31449 

0.31449 

0.31444 

0.10050 

0. 10050 

0.10046 

0.00000 

0.00000 

-0.00005 

Table  7 


TLAT  PLATE  AT  a  =  1  RAD,  CAVITATION  NUMBER  o 


Numerical 

£ 

Analytical 

K  =  20 

K  =  10 

K  =  5 

1.1 

6.324555 

6.324555 

6.324560 

6.329471 

1.3 

3.651484 

3.651484 

3.651484 

3.651730 

1.5 

2.828427 

2.828427 

2.828427 

2.828405 

1.8 

2.236068 

2.236068 

2.236068 

2.236057 

2.0 

2. 

2.000000 

2.000000 

1.999995 

Table  8 


FLAT  PLATE  AT  a  =  1  RAD,  CAVITY  VOLUME  V 


Analytical 

Numerical 

K  =  20 

K  =  10 

K  =  5 

1.1 

1.272656 

1.272656 

1.272668 

1.282604 

1.3 

1.327715 

1.327715 

1.327715 

1.328218 

1.5 

1.554475 

1.554475 

1.554475 

1.554295 

1.8 

1.975729 

1.975729 

1.975729 

1.975635 

2.0 

2.288818 

2.288818 

2.288818 

2.288768 

correspond  to  the  NACA  0010  thickness  form. 


c^:  ideal  angle  of  attack.  In  this  case  a.  =  1.54  x  —  -  _Q 

1  VJ  •  UO  /  / 

degrees . 

The  cavity  plots  for  the  three  elementary  problems  for  a  NACA  a  *  0.8 
camber  meanline  combined  with  a  NACA  0010  thickness  form  at  a  =  7°  are 
also  shown  in  Figure  21.  A  more  systematic  exposure  of  cavity  plots 
for  different  cavity  lengths  can  be  found  in  [22], 

Finally,  cavity  shapes  from  linear  and  non-linear  theory  [38]  for  a 
flat  plate  at  a  *4°  are  compared  in  Figure  22.  The  differences  are 
very  small  and  this  indicates  that  linear  theory  works  quite  well  for 
supercavitat ing  sharp-nosed  hydrofoils. 


square  root  singularities  of  the  integrands  at  n  *  0  and  n  ■  t.  The 
numerical  integrations  are  performed  thereafter  by  applying  Simpson's 
rule  with  K  uniform  intervals  in  @ . 

The  cavity  thickness  and  camber  distributions  are  computed  also 
numerically  by  applying  in  equations  (17.1)  and  (17.2)  the  algorithm 
described  by  the  formulas  (8.4)-(8.7). 

The  numerical  integrations  are  performed  for  a  flat  plate  and  an 
elliptic  foil.  The  results  are  compared  to  the  analytical  ones  given 
in  Section  18.  The  agreement  is  excellent  even  for  a  small  number  of 
intervals  in  the  numerical  integrations,  as  can  be  seen  in  Tables 
7-11. 

The  values  of  the  cavitation  number  cr  and  of  the  cavity  volume  V 

for  different  cavity  lengths  *■  are  evaluated  for  a  flat  plate  at 

a •  4° ,  a  NACA  0010  thickness  form  and  a  NACA  a  ■  0.8  camber  meanline  of 

maximum  camber  n  *  0.0679.  The  related  curves  are  shown  in 
max 

Figures  19  and  20.  From  those  curves  and  by  using  the  principle  of 

decomposition,  we  can  evaluate  the  o  vs.  I  and  V  vs.  I  curves  for  any 

combination  of  NACA  a  *  0.8  meanlines  and  four-digit  thickness  forms. 

For  example,  for  a  NACA  a  ■  0.8  meanline  with  maximum  camber  n  , 
r  max 

combined  with  a  four-digit  thickness  form  of  maximum  thickness  t 

max 

at  an  angle  of  attack  a  degrees,  we  will  have: 


a-a 


i  .  x 
+  a 


a-a 


V  =  V 


-  +  VC 


nmax  x 

Tmax 

0.0679  ° 

0.1 

X. 

Tmax 

0~0679  +  V 

U.l 

(19.2) 

(19.3) 


where: 

oa,V°:  correspond  to  the  flat  plate  at  a  *  b°. 

ac,Vc:  correspond  to  the  NACA  a  ■  0.8  camber  meanline  of  nmax 
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b.  Elliptic  Foil.  Consider  an  ellipse  of  semi-axes  (l,e)  at 
uniform  inflow  U^,  as  shown  in  Figure  18b.  According  to  the  linear 
theory,  the  streamwise  perturbation  velocity  on  the  ellipse  will  be 
constant  [39 ] : 


is 


(a) 


(b) 

Figure  18.  Special  Cases:  Flat  Plate  (a)  and  Elliptic  Foil  (b) 
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+  iFP  '  4^77~r^  '  -  *  +  - 


e+£  .  i  .  <i«2).  /  ibs_  .  lial .  i_ .  dia 

I  Z  2*  IT  J  1 1— U)  z-fu  I4tuz 


+  Jlzi  .  -2_  .  (1+22)  .  f  tL“_  .  .9.  ■(«■>,  .  _L_ 

I  t  2*  IT  *  1  t-<i»  Z— <il  1+0)* 


(17.5) 


for  z  >  t 

Analytical  solutions  for  gome  suoercavitating  hydrofoils 
»•  Flat  Plate.  For  this  case  we  have 


(18.1) 


where  a  “  angle  of  attack  as  shown  in  Figure  18a.  Then  the  equations 
(14.2),  (15.3),  (16.1)  and  (17.5),  by  using  the  integrals  given  in 
Appendix  A,  will  reduce  to  the  following  formulas,  also  given  in  [ 11 J 


and  [31]: 


o*2*  t  •  o 


q(z)  -  (P— R* z)  •  ;  z  <  t 


(18.2) 


(18.3a) 


q(z)  -  (P-R* z)  •  -  (P+Rz)  •  If—  ;  z  >  t  (18.3b) 


v(z)  -  -i>  (P-R*z)  •  ]j +  -|(P+R* z)  •  5  2  -  c  <!*.♦) 


where : 


(18.5) 


(18.6) 


Notice  that  (14.2)  expresses  o  as  an  integral  of  the  slope  of  the  lower 
hydrofoil  surface  weighted  by  terms  depending  only  on  the  cavity 
length. 


15.  Cavity  volume  V 

The  cavity  volume  V,  including  the  hydrofoil,  will  be  given  as: 

l  l 

V  =  /  h(x)  •  dx  =  -  a  f  x  •  q(x)  •  dx 

O  o  (15.1) 

Substituting  (13.6)  in  (15.1)  and  using  (A. 28),  (A. 2),  (A. 3)  and 

(14.1),  we  get: 


V  * 


o*t r*  £2 


-  a  •  /  y(£)*£*  •  d£ 

0 


(15.2) 


Finally,  by  using  (13.10),  (A. 28)  and  (A. 23),  (15.2)  becomes: 


*  u  an 

v-*2 


_n_ 

t-n 


1+n* 


j(n;t)  •  dn) 


(15.3) 


where  J(q;t)  is  given  in  (A. 23). 


16.  Cavity  source  distribution  q(z) 

Starting  with  (13.11)  and  by  using  (13.10),  (A. 28),  (A. 24),  (A. 25) 
and  (A. 26),  we  finally  get: 

q(z)  =  -0* (z)  +  (— •  - -= -  *  (  Jr2-1  -  z-  Jr2+1)  - 

2  2  •  ][T-  r2 


- -  •  (1+z2) 

7T 


6  M 

z+  U) 


■  •  dio 

(16.1a) 


for  z<t; 
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1  .  t  T(n)*dn  _  1  *  z  _  e  (x) 

?*  o  (l+n2)*(*-n)  *  1 +z2 


(13.9) 


Inversion  of  (13.9)  gives: 


Y(z)  -  -  |  •  (1+z2)  •  j  (  j  -  6*(n)] 


(1+ti2)*  (z-n) 


(13.10) 


Hotice  that  (13.10)  is  the  unique  solution  which  satisfies  also  the 
Kutta  condition  (13.4). 

Substituting  (13.8)  in  (13.6)  we  get: 


q(z)  --*+!••  Z*  (1+Z2)  •  f  — — ““ 

o  (z2-n2)* (1+n2) 


(13.11) 


The  foraulas  (13.10)  and  (13.11)  give  the  formal  solution  to  our 
problem  as  soon  as  we  determine  the  corresponding  cavitation  number 


14.  Cavitation  number  o 

The  cavitation  number  0  is  going  to  be  determined  from  the  closure 
condition  (13.5).  By  substituting  (13.6)  in  (13.5)  and  using  (A.l), 

(A. 2)  and  (A. 4)  we  get: 


-  •  7(C)  -  o 


(14.1) 


By  using  (13.10),  (A. 28),  (A. 22),  (A. 12)  and  (A. 13),  we  get  the 
following  general  expression  for  0; 


4-^Tt4 

*-(r2+l)' 


/  fa - i 

;  it-n  „ .  , 


(1+n2) 


(  \/rz+l  +  ir2-l’  n  )  •  (-  *  dy 


(14.2) 


4  2 

where:  r  ■  1  +  t  . 
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(14.3) 


(13.1) 


Then  equations  (11.7)— (11.10)  become: 

l 


_  !_  .  1  silhil  .  I  .  0  <  x  < 

0  2ir  7  £  -  X  2  * 


with  yOO  "  0 


for  1  <  x  <  i 


q(x)  .  _1_  .  X  y(0-dj  ,  e  (x)  ;  0  <  x  <  1 

"2  2.11  J  l  -  x 

0 


Y (1)  *  0 


j  q(x)  •  dx  *  0 


(13.2) 


(13.3) 


(13.4) 

(13.5) 


Using  similar  procedures  as  described  in  [lb]  and  [31],  we  derive 
the  solution  to  the  system  of  singular  integral  equations  (13.2)-(13.5) 
as  follows: 

Inverting  first  (3.2)  we  get: 

(l3-6) 


Notice  that  this  i6  the  unique  solution  which  behaves  like 
at  the  trailing  edge  of  the  cavity  (T.  Y.  Wu's  singularity  [43]). 
Substituting  (13.6)  in  (13.3)  we  finally  get: 


(a)  NACA  a  =  0.8  (n 


0.0679) 


max 


(c)  FLAT  PLATE 


. 

» 


i 


Figure  21.  The  Three  Elementary  Problems  for  NACA  a  =  0.8  combined  with  NACA  0010 
at  a  =  7°  and  l  =  1.4;  (a)  +  (b)  +  (c)  =  (d) . 
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PART  C.  CORRECTIONS  TO  THE  LINEAR  CAVITATING  HYDROFOIL  THEORY 


20.  Lighthill^s  correction  for  non-cavitating  hydrofoils 

The  linear  theory  for  non-cavitating  hydrofoils  has  already  been 
described  in  Section  1.  The  linear  solution  is  not  valid  at  a  region 
close  to  the  leading  edge  of  the  hydrofoil,  where  a  stagnation  point 
appears  and  the  perturbation  velocities  (u,  v)  are  not  small  anymore 
compared  to  the  freestream  velocity  U^, .  That  region  of  the 
non-uniformity  of  the  solution  is  of  the  order  of  the  leading  edge 
radius  p^  as  proved  by  Lighthill  in  [25]  where  also  a  uniformly  valid 
expression  for  the  surface  velocity  is  given  as  follows: 


q  =  (U  +  u-) 
ns  °°  w 


+  0(£  ) 


(20.1) 


where  e  is  our  small  parameter  and  u^  are  the  linear  streamwise 
perturbation  velocities  on  the  upper  (+)  and  the  lower  (-)  surfaces  of 
the  hydrofoil  as  shown  in  Figure  23a. 

For  example,  for  a  symmetric  hydrofoil  at  an  angle  of  attack  a, 
the  upper  surface  velocity  will  be  given  as: 


q+  *  U  (1+  uT  +  •]/ 

s  "  u  r 


(20.2) 


where  uwT  is  the  linear  streamwise  perturbation  velocity  due  to  the 
thickness  distribution  and  given  by  (1.7). 

The  minimum  pressure  point  will  be  on  the  sunction  (upper)  surface 
of  the  hydrofoil  where  the  maximum  of  qg+  occurs.  For  sufficient 
angle  of  attack,  the  minimum  pressure  point  will  be  close  to  the 


leading  edge,  at  a  chordwise  location  x  .  .  In  the  case  of  a 

min 

symmetric  hydrofoil,  x  .  will  correspond  to  the  point  which 

min 

maximizes  the  RHS  of  equation  (20.2).  For  6lowly  varying  uwT  at  the 


leading  edge  region,  which  is  the  usual  case,  x  .  is  given  as: 

mm 

xmin  ^  .  a2  (20.3) 

where  u  T  is  taken  at  x  *  0.  The  maximum  surface  velocity  q+ 

w  a ,max 

is  given  from  (20.2)  for  x  *  x  .  .  The  inception  cavitation  number 

min 

will  then  be: 

inc 

+  2 

a.  =  (q+  )  -  1 

me  s.max  (20.4) 

The  values  for  a^nc  given  by  (20.4)  agree  very  well  with  the  ones 
predicted  by  non-linear  theory  [5  ].  This  indicates  that  linear  theory 
with  the  Lighthill's  correction  gives  sufficiently  good  results  for  the 
surface  velocity  for  hydrofoils  of  practical  applications. 


21.  Lighthill^s  correction  for  cavitating  hydrofoils 

Suppose  now  that  the  hydrofoil  has  a  cavity  of  length  £  on  its 
upper  surface  as  shown  in  Figure  23b.  Then  the  leading  edge  radius  of 
the  cavity-hydrofoil  remains  the  same,  as  proven  in  Appendix  B. 
Therefore  the  velocity  on  the  cavity  surface  will  be  given  as  follows: 


,+  =  (U  +  u> 


x 


x+r 


0  <  x  <  £ 


(21.1) 


where  u£+  is  the  linear  streamwise  perturbation  velocity  on  the 
cavity  surface.  According  now  to  the  non-linear  theory,  the  velocity 
on  the  cavity  has  to  be  [38 3: 


%  =  'Z1  +  *  u»  (21.2) 

Combining  (21.1)  and  (21.2),  we  end  up  with  the  following  dynamic 
boundary  condition: 


85 


-  1]  ;  0  <  x  <  «. 


(21.3) 


I  PL 
x+  y- 


<-v  I  tfS-l-r* 

Equation  (21.3)  require*  uc+  to  be  square  root  singular  at  the 


leading  edge,  matching  asymptotically  to  a  constant  value  over  the  main 


body  of  the  cavity.  In  fact  this  is  equivalent  to  a  virtually  varying 


cavitation  number  over  the  cavity  length. 


The  kinematic  boundary  conditions  will  be  the  same  as  (2.2)  and 


(2.3).  Following  a  similar  analysis  as  described  in  Section  2,  we  end 


up  with  the  following  equations  for  the  cavity  source  and  vortex 


distributions  q(x)  and  y (x) : 


.aM  +  L  .  i  ittMC 

2  2ir  ~ 


0;  0<x<l 


(2.21) 


M  _  i_ 

2  2ir 


-  u+  -  u+ 
5-x  c  w 


0  <  x  <  l 


(21.4) 


Y(l)  -  finite 


(2.23) 


/  q(5).d£  -  0 
0 


(2.25) 


Notice  that  equations  (2.21),  (2.23)  and  (2.25)  are  the  same  as  in 


linear  theory.  Equation  (21.4)  can  be  written  as: 


U  •  |[l+o*  [l/xH  2  -  1]  -  u+ 


(21.5) 


The  left-hand  side  of  equation  (21.5)  is  in  a  similar  form  as  that  of 
equation  (2.26).  Therefore  the  coefficient  of  7  in  the  left-hand  side 


of  (21.5)  will  be  given,  according  to  (5.8),  as  follows: 


.n 


8  u+-  U  •  ^L+o*  [  I- — — - 1  ] 

(  ff*-l)«  J  w--- - -  • 

o  U 

00 

The  cavitation  number  a,  from  linear  theory  ia: 

s  ,  + 

2*u 

oL  -  /  •  f(n;s)*dn 

o  “ 

By  using  now  (7.21)  and  (21.7),  (21.6)  renders: 


f(n;s)«dn  (21.6) 


(21.7) 


a  -  [ 


a(pL;i) 


2 

]  -  1 


(21.8) 


vhere  a(p^;t)  is  defined  as: 


a (pl;A)  -  J 

o 


x+  -7 


f  (n;s)*dn 


(21.9) 


Equation  (21.8)  gives  the  cavitation  number  a  in  terms  of  the  linear 
cavitation  number  and  the  leading  edge  correction  factor 
a(pL;z).  Some  graphs  of  a(pL;i)  as  a  function  of  l  for  different 
p^  are  shown  in  Figure  24.  The  values  of  a(-p^;l)  have  been 
evaluated  by  carrying  out  the  numerical  integration  in  (21.9)  after 
making  the  transformation  (8.1).  Notice  that  as  pT  increases  with  l 
fixed  a(p^;A)  also  increases  and  therefore  a  decreases,  according  to 
(21.8).  This  implies  that  for  the  same  o>  the  bigger  the  leading  edge 
radius,  the  smaller  the  cavity  it  is,  something  which  is  also  observed 
in  the  experiments  described  in  Part  D. 

After  the  cavitation  number  o  has  been  determined  by  using 
equation  (21.8),  the  right-hand  side  of  (21.5)  is  known.  Therefore  we 
can  proceed  with  the  solution  of  the  system  of  integral  equations 
(2.21),  (21.5)  subject  to  the  restrictions  (2.23)  and  (2.25).  The 
left-hand  side  of  those  equations  is  in  the  same  form  as  in  the  linear 
theory  described  in  Section  2.  Therefore  the  solution  and  the  results 
will  be  given  by  the  same  formulas  (8.3),  (6.1)  and  (6.5)  where  instead 
of  u^+  and  <j,  the  negative  of  the  right-hand  side  of  (21.5)  and 


*  -  •  »  •  *1  e  „  » t  *■  •  j  m  «  “  — A  * *■  *  '«  *. »  *  m  '  jl  ■ «,  *  •  tm  "  »  *  »  *  %  *  • 


5©v5W9?Ss: 


2.( Vl+o-1)  must  be  put  respectively.  The  numerical  integrations  are 
performed  by  following  the  same  method  described  in  Section  8.  As  an 
application  the  effect  of  the  leading  edge  radius  on  the  cavitation  of 
a  flat  plate  is  considered.  In  this  case  u^+  is  given  by  (7.3). 

The  characteristic  curves  for  a  flat  plate  at  a  *  4  degrees  and 
different  leading  edge  radii  are  shown  in  Figure  25.  The  cavity  plots 
are  also  shown  for  a/o  *  0.06.  Notice  the  drastic  reduction  of  the 
cavity  volume  as  we  increase  the  leading  edge  radius.  This  is  of  the 
same  magnitude  as  that  predicted  by  Tulin  and  Hsu  [36]. 

Finally,  the  cavity  plots  for  three  NACA  16  thickness  forms  at  “  * 
4°  and  o  ■  1.07  are  shown  in  Figures  26  and  27,  as  predicted  by  pure 
linear  theory  and  by  the  linear  theory  including  the  non-linear  leading 
edge  corrections.  The  following  conclusions  can  be  drawn  from  the 
comparison  between  Figures  26  and  27: 

a)  Pure  linear  theory  predicts  that  by  increasing  the  thickness 
of  the  hydrofoil,  the  cavity  also  increases  for  constant  flow 
conditions  (i.e.,  a  anda  ). 

b)  Linear  theory,  including  the  leading  edge  radius  corrections, 
predicts  that  by  increasing  the  thickness  and  thus  the  leading  edge 
radius,  the  cavity  decreases  for  constant  flow  conditions.  This  is  the 
same  trend  as  that  observed  in  the  experiments  described  in  Part  D. 

22.  Effect  of  the  location  of  the  cavity  detachment  point  on  the 
cavity  shape 

Looking  in  Figure  27,  we  already  pointed  out  that  as  we  increase 


the  thickness  of  the  hydrofoil,  the  cavity  becomes  smaller  and  smaller 
for  constant  flow  conditions.  We  might  anticipate  that  as  we  go  on 
increasing  the  thickness  of  the  foil  we  will  eventually  reach  the 


on  situation  after  which  no  cavity  will  appear  on  the  foil, 
n  be  proven  by  taking  the  limit  of  equation  (21.8)  as  £-*■  0.  By 
21.7),  (21.9),  and  (7.21),  we  get: 


t(l*C  (*)>  •  j-2^  }}2  -  1  ;  as  £  -  0  (22.1) 

o  ■  [qg(0)]2  ~  1  (22.2) 

r  words,  the  inception  cavitation  number  corresponds  to  the 
>  velocity  at  the  leading  edge  for  the  non-cavitating  problem, 
rect  cavitation  number  is  given  by  (20.4)  and  corresponds  to  the 
i  surface  velocity  at  x  ■  x^. n .  Therefore  we  will  predict 
;ht  cavitation  number  if  we  start  the  cavity  at  x  *  x  .  . 
le  formulation  of  the  linear  partially  cavitating  hydrofoil 
i  with  the  cavity  starting  after  the  leading  edge  will  be  similar 
one  described  in  Section  2.  If  the  cavity  starts  at  x  *  lc  ,  the 
ry  conditions  will  become: 
xe  kinematic  boundary  conditions: 

(x,0+)  -  U  •  — ~  ;  0  <  x  <  £  and  £  <  x  <  1  (22.3) 

:  •  3x  o 

,(x,0  )  -  Vj  ;  0  <  x  <  1  (22.4) 

e  dynamic  boundary  condition: 


uc(x,°+)  -  f  •  ;  £q  <  x  <  £  (22.5) 

ity  source  and  cavity  vortex  distributions  q(x)  and  y(x)  will 
'  now  the  equations: 


v(Q.dS 

C-x 


0 ;  0<x< 1 


(22.6) 
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NL 


(22.7) 


where : 


lOO.  -  -L-  1  3.(SI-  d.$  ,  £  tj 

2  2ir  *  J  C-x  2  < 
o 

for  £  <x<£ 
o 

Y(l)  ■  finite 


-uw(x,0+) 


(22.8) 


£ 

/  q(£).d£  -  0 

£ 

o 


q(x)  “  U»-S; 


£  <x<£ 
o 


(22.9) 


(22.10) 


q(x)  -  0;  0<x<£q  and  £<x<l  (22.11) 

We  observe  again  that  the  solution  o,  y(x),  q(x)  to  the  system  of 
equations  (22.6)-(22.9)  is  unique  and  depends  only  on  the  distribution 
of  the  perturbation  velocity  uw(x,0+)  for£  <x<-£.  The  decomposition 
of  the  problem  into  elementary  problems  as  described  in  Section  3  still 
applies  as  long  as  we  keep  the  start  and  the  end  of  the  cavity  fixed 
throughout  all  those  problems.  By  using  a  similar  procedure  as  the  one 
described  in  Section  4,  we  derive  the  source  distribution  q(z)  as 
follows : 


q(z) 


1  (1+z2)  \/2~So  fS  \l s~n 

it  •  z  -Ys-z  •  *  • 


o  1+n 


.  ,  o.U 

1  /  + _ ° 

z -n  *  w  2 


)  .  di 


(22.12) 


for  s  <z<s 
o 


where : 


o  V  l-£ 


(22.13) 


The  cavitation  ruder  <?,  as  proven  in  Appendix  C,  can  be  expressed 


as: 


94 


vhcre : 


8  2*u+(n) 
f  w 

i  n 


•  f(n;s,sQ)  •  dn 


(22.14) 


-  (i+Jv  '  '  1  c53  "  c5  "ll 


(22.15) 


The  constants  C. .  C. ,  C_  are  functions  of  s  and  s  defined  in 
3’  4’  5  o 

Appendix  C.  The  function  f(n;  s,  sq)  again  satisfies  the 
restriction: 


J  f(n;s,s  )*dn  -  1 


(22.16) 


Notice  that  the  foraulas  (22.12),  (22.14)  and  (22.15)  reduce  to  (4.5), 
(5.8)  and  (5.9)  respectively  as  so*>-  0  (i.e.,  1  ♦  0) . 

The  cavity  source  distribution  can  also  be  vritten  as: 


1  1+z2 


z-s 

' — * .  / 

a— z  * 


ri-so  1+n-' 


+  [o*^  -  2-u+(z)]  •  ±  •  1/^ 


C2  Cl4lr 

tr  -  4-  * z] 


dn  + 


(22.17) 


where  the  constants  C. ,  C,  are  functions  of  s  and  s  defined  in 

Li  O 

Appendix  C. 

For  general  shape  hydrofoils,  the  expressions  (22.14)  and  (22.17) 
are  evaluated  by  making  the  transformation: 


n  -  s  +  (s-s  )  •  sin2(|-) 

O  O  i 


with  s  <,  n  £.  s  and  0  £  0  s,  ir 
o 
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(22.18) 


Then  the  involved  integrations  are  performed  numerically  by  using 

» 

Simpson's  rule  in  6 . 


The  cavity  shape  is  found  from  the  numerical  integration  of 
(22.10)  by  using  the  algorithm  described  by  the  equations  (8.4),  (8.5) 
and  (8.6),  with  $  defined  as  follows: 


5  *  *0  +  (A-*o)  •  sin2(|) 

^0<K<.H  and  0  <  $  ^  tt  (22.19) 

As  an  application,  the  partial  cavitation  for  a  NACA  16-015 

thickness  form  at  a  *  4°  has  been  considered,  as  illustrated  in 

Figure  28.  The  cavity  shapes  are  shown  for  different  starting 

positions  x  *  £q  and  the  same  end  position  x  ■  £  *  0.5.  Notice  the 

substantial  reduction  in  the  cavity  volume  and  the  cavitation  number  as 

we  move  the  start  of  the  cavity  from  x  ■  0  to  x  ■  x  as  given  by 

(20.3).  By  starting  the  cavity  farther  beyond  than  x  *  xmin»  the 

cavity  decreases  but  now  the  pressure  minimum  will  occur  on  the  wetted 

part  of  the  foil  in  front  of  the  cavity,  something  which  is  physically 

unacceptable.  In  Figure  29,  the  cavity  shapes  are  shown  for  a  NACA 

16-006  thickness  form  at  <*  *  4°  and  fort  ■  0  ,  £  ■  x  .  . 

Comparing  Figure  29  to  Figure  28,  we  observe  that  starting  the  cavity 

at  l  *  x  .  has  smaller  effect  on  the  solution,  the  thinner  the 
o  min 

hydrofoil. 

The  leading  edge  corrections  can  also  be  included  when  the  cavity 
starts  at  x  *  >0.  In  that  case,  the  dynamic  boundary  condition  on 

the  cavity  surface  will  become: 


u+  -  v  t 


x+ 


- 1) 


£  <  x  <  £ 
0 


(22.20) 


The  cavitation  number  0  will  again  be  given  by  (21.8)  where  now 


corresponds  to  the  linear  theory  with  the  same  ^ q  and  £ ,  and  where 
instead  o£  a(pT;0,  S(p_  )  is  used,  defined  as: 

L  L  O 


a(pL;«.,*0)  “  / 

s 


2  - 

— —  •  f(n;s,so)»dn 


(22.21) 


After  a  has  been  determined,  we  proceed  by  solving  the  linear  problem 

with  the  same  z  and  i,  but  with  u^+  replaced  by  the  negative  of 

the  right-hand  side  of  (21.5).  In  Figure  30,  the  cavity  shapes  for  a 

NACA  16-006  thickness  form  at  a  *  4°  and  Z  =  0.5,  are  shown  for  Z 

*  0  and  *  Xmin*  aS  Pre<*icte<*  by  li-near  theory,  including  the 

leading  edge  corrections.  Notice  the  negligible  differences  between 

the  two  cavity  shapes  and  their  cavitation  numbers.  This  means  that  in 

the  linear  theory  with  leading  edge  corrections,  the  point  where  the 

cavity  starts  has  little  effect  on  the  solution,  something  which  has 

also  been  concluded  from  the  non-linear  theory  [38],  However,  by 

letting  Z  *  x  .  ,  we  predict  the  right  cavitation 
o  min 

inception  number. 

By  comparing  Figure  30  with  Figure  14,  we  also  observe  that 
including  the  leading  edge  corrections  makes  the  cavity  thinner  and,  in 
fact,  closer  to  the  non-linear  result.  However,  the  predicted 
cavitation  number  is  still  higher  than  the  non-linear  one.  The 
predicted  cavity  shape  for  the  same  cavitation  number  as  the  non-linear 
one  is  shown  in  Figure  31a.  It  seems  that  we  still  overpredict  the 
cavity  volume  for  fixed  flow  conditions,  but  not  as  much  as  with  pure 
linear  theory.  Nevertheless,  including  the  leading  edge  corrections 
makes  the  cavity  shape  very  close  to  the  non-linear  one  at  its 
beginning.  The  discrepancies  at  the  end  of  the  cavity  can  partly  be 
attributed  to  the  different  cavity  termination  models  in  the  two 


theories.  This  is  encouraged  by  the  fact  that  non-linear  theory  using 
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a.  Cavity  plots  for  NACA  16-006  thickness  form  at  a  =  4°,  a  = 

—  Non-linear  theory  [38] 

-  Pure  linear  theory 

— -  Linear  theory  including  leading  edge  corrections 


Figure  31b.  Cavity  plots  for  flat  plate  at  a  *  5°,  a  =  1.246 

-  ■  "  Non-linear  theory — Riabouchinsky  model  [33 1 

-  Non-linear  theory — re-entrant  jet  model  [23 1 


Riabouchinsky  end  plate  model  [33]  underpredicts  the  extent  and  volume 
of  the  cavity  for  a  flat  plate,  compared  to  the  non-linear  theory  using 
the  re-entrant  jet  model  [23],  as  shown  in  Figure  31b. 


23.  Open  cavity  model 

The  open  wake  models  have  been  developed  in  order  to  simulate 
better  the  actual  flow  behind  the  cavity  and  the  ultimate  wake  defect, 
which  is  directly  associated  to  the  generated  drag.  A  review  of  these 
models  has  been  given  by  Tulin  in  [35]. 

The  linearized  version  of  the  open  wake  models  was  first 
introduced  by  Fabula  in  [7  ].  According  to  that  model,  the  wake  has  a 
constant  thickness  £  equal  to  the  thickness  of  the  cavity  at  its  end, 
as  shown  in  Figure  32a.  The  formulation  of  the  linearized  problem  will 
be  exactly  the  same  except  for  the  cavity  closure  condition  which  will 
become: 


6  =  / 
0 


q(Q 

u 


d£ 


or  by  using  the  transformation  (A. 3) 


(23.1) 


_  J  q(z)  #  2* z 


0  U 


(1+z2)2 


dz 


(23.2) 


The  cavity  source  q(z)will  still  be  given  by  (4.5),  since  no  use  of  the 
closure  condition  has  been  made  to  derive  it.  If  the  change  in 
cavitation  number  due  to  the  finite  thickness  6  at  the  end  of  the 
cavity  is  Ac  ,  then  (23.2)  with  (4.5),  (5.3),  (5.4)  and  (7.21)  will 
render: 


Ao 


(23.3) 


where  D  depends  only  on  the  cavity  length  £.  Notice  that  Aa  is 


proportional  to  6  for  a  given  cavity  length  and  also  that: 

Aa  <  0  (23.4) 

The  algorithm  to  find  the  cavity  shape  will  remain  exactly  the  same  as 
in  linear  theory,  after0  has  been  adjusted  by  (23.3). 

Including  the  leading  edge  radius  corrections,  in  the  case  an  open 
cavity  model  is  used,  is  straightforward.  The  resulting  formulas  will 
be  the  same  as  described  in  Section  21,  except  that  must  be 
changed  by  an  amount  given  by  (23.3).  If  the  cavity  starts  after  the 
leading  edge  of  the  hydrofoil,  then  by  using  the  formulas  developed  in 
Appendix  C,  we  will  get: 


The  effect  of  the  open  wake  thickness  6  on  the  cavity  shape  for  a 
NACA  16-006  at  a  =4°  is  shown  in  Figure  32b.  By  increasing  6,  the 
following  conclusions  are  drawn: 

a)  The  cavity  becomes  thicker  at  its  end  and  slightly  thinner  at 
its  forebody. 

b)  The  cavity  volume  increases  slightly  for  the  same  cavity 
length. 

c)  The  cavitation  number  decreases  substantially. 

d)  The  cavity  length  and  volume  decrease  substantially  for  the 
same  o • 

The  present  theory  with  6=0,  is  applied  for  some  NACA  16  series 
Section  forms  and  the  results  are  shown  in  Appendix  D,  in  comparison 
with  other  theories. 
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PART  D.  CAVITATION  EXPERIMENTS  OH  ROUND-NOSED  HYDROFOILS 
14.  Purpose  and  description  of  the  experiment 

In  section  21  ve  concluded  that  the  leading  edge  radius  plays  a 
lignificant  role  in  the  partial  cavitation  of  a  hydrofoil.  Namely,  ve 
found  that  the  cavitation  number?^  resulting  from  pure  linear  theory 
could  be  corrected  due  to  the  leading  edge  radius,  as  follows: 


[ 


a(pL;i) 


2 

]  -  1 


(21.8) 

It  would  be  interesting  to  see  how  well  this  predicts  what  happens  in 
reality.  The  MIT  Hydrodynamics  Laboratory  has  at  its  disposal  three 
round-nosed  aluminum  hydrofoils  with  identical  projected  area,  camber 
distribution,  and  maximum  thickness,  but  with  different  leading  edge 
radii.  Therefore,  by  testing  those  hydrofoils,  ve  would  be  able  to 
study  solely  the  effect  of  their  leading  edge  radii  on  their 
cavitation. 

The  three  hydrofoils  have  the  following  main  characteristics  (with 
respect  to  the  chord  length): 


No.  1:  This  is  a  NACA  a  *  0.8  meanline  of  n 


max 


■  0.02  combined 


with  a  VLR  thickness  form  [20]  of  t _  ■  0.04,  p 


(1) 


4  x  p 


(2) 


0.002867  and  x, 


L  .  .  T.E. 

No.  2:  This  is  a  NACA  a  ■  0.8  meanline  of  n 


max 

-  0.004. 


max 


with  a  NACA  66  thickness  form  of  t  •  0.04,  p_ 

max  ’  L 


0.02  combined 

(2) 


•  0.000717  and  t_  _  -  0.004. 

T  •  £  • 

No.  3:  This  is  a  NACA  a  ■  0.8  meanline  of  n 
with  a  VLR  thickness  form  of  t 


2.25  x  p 


(2) 


max 


max 

0.04,  p 


0.02  combined 
(3) 


-  0.001613  and  x_  „  -  0.004. 
T.E. 


Finally  it  was  decided  to  test  the  hydrofoils  No.  2  and  No.  3, 
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(n)  The  reduction  in  the  cavity  size  by  increasing  the  leading 
edge  radius  of  the  hydrofoil  has  also  been  confirmed  from  the 
experiments  conducted  at  the  MIT  water  tunnel. 

RECOMMENDATIONS 

The  conclusions  of  the  thesis  suggest  the  following  directions  for 
further  research: 

a)  To  investigate  the  effect  of  the  leading  edge  radius  on  the 
load  distribution  and  the  hydrodynamic  coefficients  of  a  cavitating 
hydrofoil. 

b)  To  apply  the  leading  edge  corrections  to  supercavitating 
round-nosed  hydrofoils.  In  this  case,  the  corrections  are  expected  to 
have  less  impact  on  pure  linear  theory. 

c)  To  conduct  more  experiments  on  round-nosed  hydrofoils, 
preferably  of  thick  leading  edge  radii,  and  to  compare  the  measured 
values  of  the  cavitation  numbers  for  different  cavity  lengths  with  the 
ones  predicted  by  the  proposed  theory.  When  applying  the  theory,  the 
wall  tunnel  effects  should  be  properly  taken  into  account.  The  values 
of  an  equivalent  open  cavity  thickness  that  would  make  the  present 
theory  fit  the  experiment  should  be  found,  and  their  correlation  to 
the  main  variables  of  the  problem  should  be  investigated. 

d)  To  apply  the  leading  edge  corrections  to  the  cavitating 
propeller  problem  by  incorporating  a  local  Lighthill's  correction  at 
each  spanwise  position  on  the  leading  edge  of  the  blade.  This  is 
expected  to  have  an  insignificant  effect  on  the  computing  time  of  the 
existing  code. 


super cavit at ing  hydrofoils  at  noderate  angles  of  attack. 

(g)  Linear  theory  has  been  shown  to  overpredict  substantially  the 
extent  and  volune  of  the  cavity  for  partially  cavitating  round-nosed 
hydrofoils. 

(i)  The  leading  edge  radius  of  the  hydrofoil  has  been  introduced 
by  incorporating  Lighthill's  correction  in  the  linearized  formulation 
of  the  problem.  This  led  to  the  introduction  of  a  virtual  cavitation 
number  varying  over  the  cavity  length  and,  in  particular,  square  root 
singular  at  the  leading  edge  of  the  hydrofoil.  The  implementation  of 
the  leading  edge  radius  corrections  into  our  computer  code  for 
partially  cavitating  hydrofoils  was  trivial  and  had  no  significant 
effect  on  the  computing  time. 

(j)  The  effect  of  the  leading  edge  radius  on  the  partial 
cavitation  of  a  hydrofoil  has  been  found  to  be  significant.  In 
particular,  it  was  shown  that  by  increasing  the  leading  edge  radius  of 
a  hydrofoil,  the  extent  and  volume  of  the  cavity  decreased,  for  fixed 
flow  conditions. 

(k)  The  effect  of  the  location  of  the  cavity  detachment  point, 
when  the  leading  edge  radius  corrections  were  included,  has  been  found 
to  be  negligible  for  long  cavities.  However,  by  starting  the  cavity  at 
the  minimum  pressure  point  we  predicted  the  right  cavitation  inception 
number. 

(l)  The  effect  of  the  cavity  termination  model  on  the  cavity 
extent  and  volume  has  been  found  to  be  significant. 

Cm)  The  present  theory,  compared  with  the  non-linear  theory  by 
Uhlman,  has  been  found  to  overpredict  the  effect  of  leading  edge  radius 
on  the  cavity  size,  but  not  as  much  as  the  short  cavity  theory  by  Tulin 


C01CLUSI0NS 


The  following  have  been  accomplished  during  the  evolution  of 
this  thesis: 

(a)  The  linearized  theory  for  the  partial  and  the  supercavitating 
hydrofoil  problem  has  been  stated  in  a  uniform  way  in  terms  of  integral 
equations  of  unknown  source  and  vortex  distributions  on  the  foil  and  on 
the  cavity. 

(b)  The  general  problem  for  partial  and  supercavitating  hydrofoils 
has  been  decomposed  into  three  elementary  problems,  one  of  camber,  one 
of  thickness  and  one  of  angle  of  attack.  This  has  enabled  us  to  study 
separately  the  effect  of  camber  and  thickness  on  the  cavitation  of 
hydrofoils. 

(c)  The  general  solution  of  the  linear  partially  cavitating 
hydrofoil  problem  for  a  given  cavity  length  has  been  shown  to  depend 
only  on  the  distribution  over  that  length  of  the  non-cavitating 
streamwise  perturbation  velocity  u^+  on  the  upper  hydrofoil 
surface.  Expressions  for  the  cavitation  number,  the  cavity  volume  and 

the  cavity  source  distribution  have  been  given  in  terms  of  integrals  of 

+ 

u 

w 

(d)  For  the  supercavitating  hydrofoil  problem,  expressions  for  the 
cavitation  number,  the  cavity  volume  and  the  cavity  source  distribution 
are  given  in  terms  of  integrals  of  the  slope  of  the  wetted  hydrofoil 
surface. 

(e)  Very  accurate  and  efficient  computer  codes  have  been  developed 
for  the  application  of  linear  theory  in  partial  and  supercavitating 
hydrofoils. 

(f)  Linear  theory  has  been  shown  to  be  accurate  for  sharp-nosed 


Measured 

a 

^1+0 

%<3> 

Actual 

a 

0. 

2.45 

0.1 

1.14 

0.2 

0.96 

0.3 

0.86 

0.4 

0.83 

Table  12.  Experimental  results  for  foil  No.  3 


Measured 

a 

^1+a 

%<2> 

Actual 

o 

0. 

2.944 

0.1 

1.262 

0.2 

1.047 

0.3 

0.956 

0.4 

0.91 

1.045 


0.905 


0.795 


Table  13.  Experimental  results  for  foil  No.  2 


actual  leading  edge  radii  or  by  testing  hydrofoil  No.  3  which  has  a 
larger  leading  edge  radius  and  therefore  is  less  vulnerable  to 
machining  imperfections. 
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changed  from  3.  degrees  tc*2.7  degrees,  due  to  the  angular  deflection 
of  the  hydrofoils  around  their  hinges,  produced  by  the  developed 
loading  during  the  experiments.  The  reduction  in  the  cavity  extent  due 
to  the  increase  in  the  leading  edge  radius,  for  fixed  a  /a  ,  is 
significant  especially  for  short  cavities.  A  similar  reduction  in 
cavity  volume  is  expected  but  the  actual  curves  are  not  shown,  since  it 

is  difficult  to  measure  the  cavity  volume. 

(2)  (3) 

If  a  and  o  are  the  cavitation  numbers  predicted  by  the  present 
theory  for  the  two  hydrofoils,  then  at  a  given  cavity  length  the 
following  relationship  will  hold: 


(25.1) 


(3)  (2)  (3) 

where  a  and  a  are  the  correction  factors  corresponding  to p ^  and 

(2) 

p  L  respectively,  at  the  given  cavity  length.  Equation  (25.1)  was 

(2)  (3) 

derived  from  (21.8)  under  the  assumption  that  o ^  ■ a ^  which  is 

true  since  the  two  hydrofoils  are  identical  except  for  their  leading 

edge  radii,  which  do  not  affect  the  linear  solution  anyway.  In  order 

(3) 

to  compare  the  present  theory  with  the  experiments  we  evaluated  o 

(2  ) 

from  (25.1)  by  using  the  experimental  values  of  aK  .  The  respective 
l  VB.a/a  curve  is  also  shown  in  Figure  34.  We  observe  that  the 
present  theory  overpredicts  the  reduction  in  the  cavity  extent,  for 
fixed  a /a  .  This  may  be  attributed  to  any  differences  of  the  actual 


leading  edge  radii  from  the  designed  ones  due  to  machining 
imperfections,  especially  for  hydrofoil  No.  2  which  has  the  smaller 
leading  edge  radius.  This  could  be  verified  by  either  measuring  the 


The  cavitation  number  for  each  cavity  length  is  evaluated  by 
substituting  in  formula  (2.1)  the  recorded  values  for  ^  and  . 

The  resulting  numbers  for  the  tvo  tested  hydrofoils ,  designated  as 
measured  cr  are  shown  in  Tables  12  and  13.  The  corresponding  values 
at  l  m  0  are  the  cavitation  inception  numbers. 

We  have  also  been  able  to  measure  the  flovfield  velocity  on  the 
cavity  surface,  using  the  laser  Doppler  anemometer  installed  at  the  MIT 
water  tunnel,  operating  in  forward  scatter  mode.  In  fact,  by 
approaching  the  cavity  with  the  laser  beams  from  above,  we  recorded  the 
measured  velocity  just  before  the  signal  was  lost  by  the  tracker, 
indicating  that  the  beams  entered  the  cavity.  The  measured  cavity 
velocities  were  found  to  be  independent  of  the  spanwise  or  chordwise 
location  on  the  cavity.  Their  values  relative  to  the  freestream 
velocity  are  shown  at  the  fourth  column  of  Tables  12  and  13.  We 
observe  that  they  are  smaller  than  the  values  of  |1  ♦  o  based  on  the 
measured  a  ,  also  shown  at  the  third  column  of  the  same  Tables.  The 
values  of  the  cavitation  numbers  which  would  predict  the  same  cavity 
velocities  as  the  measured  ones,  designated  as  actual  cr's,  are  also 
shown.  The  difference  of  the  measured  from  the  actual  cavitation 
number  varies  slightly  with  the  cavity  length  and  finally  is  attributed 
either  to  the  entrapped  air  inside  the  cavity,  or  to  calibration  errors 
of  the  manometers.  For  any  comparisons  between  experiments  and  theory 
the  actual  cavitation  numbers  should  be  used. 

In  Figure  34  the  i  vs.  a  /a  curves  are  shown  for  the  two 
experiments.  The  actual  angle  of  attack  was  measured  and  was  found 
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The  hydrofoils  were  put  in  the  twenty  by  twenty  inches  test 
section  of  the  MIT  variable  pressure  water  tunnel  [18],  as  shown 
schematically  in  Figure  33.  They  were  mounted  from  their  hinges  H  on 
specially  provided  holes  on  the  vertical  windows  of  the  water  tunnel. 
The  angle  of  attack  was  set  equal  to  3  degrees  for  both  tests,  by 
rotating  and  then  locking  the  hinges  via  devices  mounted  on  the  outer 
surface  of  the  windows.  Thin  rubber  layers  were  placed  in  between  the 
tips  of  the  hydrofoils  and  the  inner  surfaces  of  the  vertical  windows 
to  prevent  any  cross  flow  due  to  the  pressure  differences  between  the 
lower  and  the  upper  surfaces  of  the  hydrofoil.  Marker  lines  were  drawn 
on  the  upper  surface  of  the  hydrofoil  in  order  to  facilitate  the  cavity 
length  readings. 

During  both  tests  the  flow  velocity  has  been  kept  constant  around 

19.5  ft/sec.  By  changing  p  ,  we  changed  the  cavitation  number  a  and 

00 

therefore  the  size  of  the  cavity.  The  following  observations  were 
made: 

a)  The  cavities  appeared  very  two-dimensional. 

b)  The  end  of  the  cavity  was  oscillating  around  some  mean 
position.  The  longer  the  cavity,  the  bigger  the  oscillations. 

Cavities  longer  than  60%  of  the  chord  fell  quickly  into  the  3/4  chord 
instability  regime  in  which  the  cavity  was  moving  violently  back  and 
forth. 

c)  The  cavity  surface  at  its  beginning  was  smooth  and 
transparent,  becoming  rough  and  opaque  towards  the  end,  followed  by  a 
wake  of  similar  structure. 


[43]  Wu,  T.  Y. ,  "A  Simple  Method  for  Calculating  the  Drag  in  the 
Linear  Theory  of  Cavity  Flows,"  California  Institute  of  Technology, 
Report  Ho.  85-5,  August  1957. 


APPENDIX  A 


List  of  Integrals 

In  this  appendix,  a  list  of  integrals  used  throughout  our  analysis 
is  given,  along  with  some  instructions  for  their  derivation. 

The  first  integrals  (A.1)-(A.7)  have  been  evaluated  by  using  the 
transf ormat ion : 

2  0  2  0  * 
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The  next  integrals  (A.8)-(A.18)  have  been  taken  from  [31]  where  they 


have  been  evaluated  by  straightforward  application  of  the  calculus  of 
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Whenever  ve  have  to  change  the  order  of  integration  in  a  double 
singular  integral,  ve  use  the  Poincare— Bertrand  formula  l27l: 
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APPENDIX  B 


Leading  edge  radius  for  cavitatine  hydrofoils 
The  point  where  the  cavity  starts  for  the  case  of  round-nosed 
hydrofoils  can  be  determined  in  non-linear  theory  by  the 
Villat-Brillouin  condition  [4,  AO],  according  to  which  the  curvature  of 
the  cavity  and  foil  should  be  continuous  at  that  point.  Some 
semi-empirical  criteria  to  take  into  account  viscous  effects  can  also 
be  used  [2  ].  For  the  case  of  sharp-nosed  hydrofoils,  the  cavity  will 
start  from  the  leading  edge  tangent  to  the  lower  hydrofoil  surface, 
going  ahead  of  the  foil  for  a  short  distance  and  then  turning  back 
towards  the  trailing  edge  of  the  hydrofoil. 

In  the  linear  theory  described  in  Parts  A  and  B  of  this  thesis,  we 
assume  that  the  cavity  starts  at  the  leading  edge  of  the  hydrofoil  for 
either  sharp-  or  round-nosed  hydrofoils.  In  both  cases  the  asymptotic 
behavior  of  q(x)  as  x*-0  can  be  found  by  using  (4.5)  and  (4.3)  as: 


,  x  .  (const) 
q  (x)  'v  - — ; -  as  x-n) 

Xi 


(B-l) 


which  by  using  (2.19)  will  give: 


dh  _  (const) 

~ s — as  <«-2> 

For  round-nosed  hydrofoils  the  slope  of  the  upper  surface  behaves 
like: 


dx  V  2  ' 


as  x-O 


where  is  the  leading  edge  radius  of  the  hydrofoil, 
surface  nc(x)  can  be  expressed  as: 


(B-3) 

But  the  cavity 
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and  the  slope  of  the  cavity  surface  will  be  given  as  the  sum  of  the 
foil  and  the  cavity  thickness  slopes: 


The  asymptotic  behavior  of 


dn  dri 

_ c  m  _ u 

dx  dx 


3ric 

-r —  will  be: 
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dh 

dx 


(B-5) 


since  the  second  term  in  the  RHS  of  (B-5)  is  subdominant  according  to 
(B-2)  and  (B-3).  Therefore  the  slope  of  the  cavity  is  the  same  as  the 
slope  of  the  hydrofoil  at  the  point  where  the  cavity  starts.  Similarly 
we  may  conclude  that  the  radius  of  the  curvature  of  the  cavity  surface 
at  its  start  is  p^»  the  same  as  the  curvature  of  the  hydrofoil  at 
that  point. 

For  supercavitating  round-nosed  hydrofoils,  if  the  cavity  starts 
at  the  leading  edge  of  the  hydrofoil,  then  the  formula  (16.1a)  will 
render  [22]: 


_  * 

q(z)  'v*  -2.9  (z)  as  z -+0  (B-7) 

which  also  means  that  the  cavity  has  the  same  leading  edge  radius  as 
the  hydrofoil. 

When  the  leading  edge  corrections  are  included,  again  the  leading 
edge  radius  remains  the  same  as  that  of  the  foil,  since  the  cavity 
source  distribution  will  have  again  a  fourth  root  singularity  at  the 
leading  edge,  as  can  be  concluded  from  Section  21. 


APPENDIX  C 


Formulas  used  in  Section  22 

The  cavity  source  distribution  in  the  case  the  cavity  starts  at  £o>0 
has  been  given  in  Section  22  as  follows: 
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where : 
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The  cavitation  number  a  is  going  to  be  determined  from  the  cavity 
closure  condition  (22.9),  which  in  terms  of  z  will  become: 
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Substituting  (22.12)  in  ( C— 1 )  and  using  (A. 28)  we  get: 
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The  integral  I(n;  s,  sq)  is  identical  to  I(T1;  s)  defined  in  (5.7), 
for  ■  0.  It  can  also  be  expressed  as: 
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where  C^,  are  constants  depending  on  s,  sq  and  defined  as: 
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The  cavity  source  distribution  q(z)  can  also  be  put  in  the  following 
form,  starting  from  (22.12): 
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The  constants  Cj,  Cj,  C^,  for  so  ■  0  reduce  to  the 

constants  defined  by  the  expressions  (A. 15),  (A. 14),  (A. 9),  (A. 8). 

Similarly  Cj  -*,D.x  as  sq*^  0.  Analytical  expressions  for  , 

Cj ,  Cj,  C4,  Cj  as  functions  of  s  and  sq  can  be  found  through 
a  cumbersome  application  of  the  calculus  of  residues.  Instead  we 


perform  the  integrations  involved  in  their  definitions  numerically 
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after  making  the  transformation  (22.18) 


APPENDIX  D 


Characteristic  curves  and  cavity  plots  for  NACA  16  thickness  forms 

The  computer  codes  developed  for  the  application  of  the  theory 
described  in  Parts  A  and  C  of  the  thesis  have  been  run  for  symmetrical 
NACA  16  sections  with  thickness/chord  ratios  of  6%,  9%,  and  12%,  at 
an  angle  of  attack  a  ■  4°.  The  input  £  and  P^  values,  in  the 
case  the  leading  edge  corrections  are  included,  are  shown  in  Table  Dl . 
The  values  for  £  correspond  to  x  .  and  they  have  been  evaluated 
by  using  (20.3).  The  results  have  been  found  to  be  as  accurate  and 
convergent  as  stated  in  Section  8.  The  required  CPU  time  has  been 
negligible. 

The  characteristic  curves  as  given  by  pure  linear  theory,  linear 
theory  including  leading  edge  corrections,  non-linear  theory  [38]  t  and 
short  cavity  theory  [36]  are  shown  in  Figures  D1-D10.  The  following 
conclusions  can  be  drawn: 

a)  Pure  linear  theory  predicts  that  as  the  foil  thickness 
increases,  the  cavity  length  and  volume  also  increase  for  fixed  flow 
condit ions. 

b)  Linear  theory,  including  the  leading  edge  corrections, 
non-linear  theory  and  short  cavity  theory,  all  predict  that  by 
increasing  the  foil  thickness,  the  cavity  length  and  volume  decrease 
for  fixed  flow  conditions.  The  reduction  in  cavity  volume  as 
predicted  by  the  present  theory  is  higher  than  that  predicted  by  Uhlman 
[38],  but  not  as  high  as  that  by  Tulin  and  Hsu  [36], 

c)  Comparing  Figure  D5  with  D6,  we  observe  that  for  fixed  cavity 
length,  the  present  theory  predicts  a  slightly  higher  cavity  volume 
than  the  non-linear  theory.  Comparing  also  Figures  D2  with  D3,  we 
conclude  that  it  is  the  difference  between  the  corresponding  cavitation 
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numbers  which  makes  the  two  theories  different.  By  adjusting  the 
thickness  of  the  cavity  at  its  trailing  edge,  as  described  in  Section 
23,  we  could  very  well  make  the  two  theories  produce  equivalent 
results.  However,  since  the  cavity  termination  model  in  non-linear 
theory  is  not  necessarily  physically  correct,  it  is  still  necessary  to 
make  use  of  experimental  results  to  judge  which  theory  is  better. 

Finally,  some  cavity  plots  for  the  three  NACA  sections,  as 
predicted  by  linear  theory  including  leading  edge  corrections,  are 
shown  in  Figures  D11-D13. 


Table  Dl.  Input  data  for  the  NACA  16  thickness 
forms  at  a  =  4° 


1. 


Figure  D2.  i  vs.  a/a  for  NACA  16  series  sections  of  6%,  9%,  and 

12%  thickness/chord  ratios  at  a  =  4°  using  linear  theory 
including  leading  edge  corrections. 


£  vs.  a/a  for  NACA  16  series  sections  of  6%,  9%,  and 
thickness/chord  ratios  at  a  =  4°  using  non-linear 
theory  [38]. 


Figure  D5.  V/a  vs.  £  for  NACA  16  series  sections  of  6%,  9%,  and  12% 
thickness /chord  ratios  at  a  -  4°  using  linear  theory 
including  leading  edge  corrections. 


Figure  D7.  V  vs.  a/a  for  NACA  16  series  sections  of  6%,  9%,  and  12% 
thickness/chord  ratios  at  a  •  4°  using  pure  linear  theory. 
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V  vs.  a/o  for  NACA  16  series  sections  of  6%,  9%  and  12% 
thickness/chord  ratios  at  a  *  4°  using  non-linear  theory 
[38]. 
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Figure  DIO.  V  vs.  a/a  for  NACA  16  series  sections  of  6%,  9%,  and  12% 
thickness/chord  ratios  at  a  =  4°  using  short  cavity 
theory  [36]. 
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